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INTilODUCTIOM. 


rilUIS Essay aims at a complete c*alciilatioii of the effects pro(hiC(‘(l by the 
J- action of a planet on the motion of the moon under the following 
limitations and conditions: 

(1) The problem of the motion of the moon under the action of the 
sun (supposed to move I’ound tht‘ ceiitrt* of mass of tlu‘ earth and moon in 
a fixed t‘lliptic orbit) and the t‘arth, is considered to have been completely 
solved. 

(2) All the bodi(\s are supposed to attract in tht‘ same manner as 
particles of masses <*(pial to their actual masses and situatc‘d at the centres 
of mass. 

(3) All the planets are supposed to moves in fixed elliptic orbits, 
i.e.f the effect of the action of a planet transmitted either through the earth 
or through another planet is neglected. 

(4) Perturbations of the first order with respect to the ratio of the 
mass «)f a planet to that of the sun are alone calculated. 

(5) The exception to the above limitations occurs in the periods of 
revolution of the apse and node of the moon’s orbit. The.se periods are n<jt 
exactly those arising from (1) but they are the observed periods or, what 
amounts to the same thing owing to the clo.se agreement between the 
observed and calculated periods, the pericxls after all known causes have been 
included. The point is only of importance in terms of very long period. 

(G) All coefficients greater than (/'’Ol in longitude, latitude and 
parallax have been obtained. Many are also given which are le.ss than 
0" 01 whenever they have been accurately calculated. There are, in addition, 
classes of terms of short period which run in series and which in the 
aggregate will add up at certain times to much more than 0 *01 : these have 
been found to be 0'’ 002. 
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(7) The maximum period considered is 3500 years, but as the sieve 
in Section iv. retained a few terms of longer period, these were also included 
in the general scheme. 

The methods here adopted have been constructed mainly to overcome 
the difficulties which have in the past prevented an accurate computation 
of the long period terms. They were, however, found to be equally useful 
for finding the terms which have periods of a year or less. These difficulties 
include the development of the parts of the disturbing function which 
depend on the coordinates of the earth and planet; the accurate calculation 
of the derivatives of the moon’s coordinates with respect to n, the moon’s 
mean motion ; the uncertainty arising from the possible omission of terms 
of long period ; and the frequent appearance of small coefficients as the 
difference of two large numbers. 

In Section i., the equations of variations for the lunar elements have 
been recomputed with the use of a semi-canonical system of elements. The 
equations for the ordinary system of elements were first given by G. W. Hill 
and independently, though at a later date, by S. Newcomb ; they were 
recalculated in Hill s form by R. Radau. In the present system, errors 
arising from the slow convergence of Delaunay’s series have been avoided ; 
in fact his literal expressions have only been used in small terms where 
derivatives with respect to n were required but where the maximum possible 
error could make no difference in the final results. 

In Section ii. Hill’s method of dividing the disturbing function into 
a sum of products in which the first factor of each product is independent 
of the lunar cooi*dinates and the other of the planet’s coordinates, is exhibited 
as part of a general theorem. 

By referring these coordinates to the ti^ue place of the sun’s radius vector, 
I have obtained the first factors directly from the expansion of the inverse 
first power of the distance between the planet and the earth (1/A) and of its 
derivatives with respect to certain of the elements of the earth and planet. 
Only one expansion is therefore required for all these factors, namely, that 
of 1/A,and this has been given by Leverrier in a literal form in powers of the 
eccentricities and mutual inclination. The expansion also contains the 
coefficients in the expansion of I/Aq (the value of a^A when the eccentricities 
and inclination are zero) and its derivatives with respect to the planet’s mean 
distance : the formulae for finding the coefficients and their derivatives are 
here put into forms which admit of rapid and simple computation. 

The factoi's containing the moon’s coordinates, together with their deriva- 
tives with respect to all the lunar elements except n, are found from the 
results of my lunar theory; a special method which I gave five years ago and 
which does not require the use of literal series in powers of n'jn has been 
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used to find the derivatives with respect to n. These methods for finding 
the planetary and lunar factors are set forth in Section ill. 

The only one of the difficulties previously mentioned which has not been 
considered up to this point is the danger of omitting long-period terms with 
sensible coefficients. In Section iv. formulae are constructed which permit 
one to find rapidly an upper limit to the magnitude of any coefficient. By 
means of them, all terms having coefficients greater than 0"'01 and periods 
less than 3500 years have been sifted out; there are about 100 of such terms, 
excluding the terms of short period for which no sieve was required. 

Sections v., vi. consist of numerical results. It may be noted that the 
values of Apy Bp, ... are also required for finding the perturbations of the 
earth by the planets and of the planets by the earth, while those of Mi and 
of their derivatives (as well as the equations of variations contained in 
Section I.) are available for the computation of lunar perturbations other than 
those due to the direct action of the planets. 

No new inequalities sufficiently great to account for the observed dis- 
crepancy between theory and observation appear from the direct action of 
the planets, as shown in the tables of Section vi. Radau^s well-known list 
of terms in longitude has required considerable extension as far as the short 
period terms are concerned, and a few new long period inequalities with 
small coefficients have been computed. The more extensive developments 
of this essay have shown that some of his coefficients re(iuire alteration, but 
there is a general agreement for all those portions which he has taken into 
account. 


Only a few slight verbal changes and corrections of errors in copying have 
been made to the first five sections, with two exceptions mentioned below, 
since the award of the examiners. But I have gone over all the computations 
for finding the short period terms and the larger long period terms during 
the year that has elapsed and have made the following corrections to the 
results in Section vi. : 


Argument 

i+32’-10r+33'‘ 
-^■-16^+ 187-151” 
29^-267+ 112” 
2D-/ + 21^-207-87” 


Former coefticient 
- 0"-35 
-15"-22 
4- 0"-117 
+ 0"-lll 
+ 0"-040 


Corrected coefficient 
+ 0"*36 
-14" *55 
+ 0"-108 
+ 0"-126 
- 0" 038 


together with their accompanying short period secondaries. 

The signs of those coefficients containing h on p. 86 have been changed. 
The annual mean motion of the perigee has been altered from 2"*66 
to 2"*69. 

* I gave this value in a paper referred to on page 3 below. 
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The values of and of its derivatives with respect to ti, e, k under the 
argument 0 on page 61 have required a factor 2. This error necessitated 
slight changes in some of the coefficients whose primaries were independent 
of the lunar angles; the largest correction was one of — 0"‘019 in the coeffi- 
cient of sin (T - F). 

A wrong sign in the computation of the equations of variations (see 
Errata at the end of the volume) gave rise to a few almost insensible 
changes in certain coefficients. 

The additions are : 


Argument 

Coefficient 

4M-2T+6r 

-0"012 

21-^21) + 4 

4-0"‘017 

2/-2D + 8J/-6f+63“ 

-0"019 

2/-2D + 83f-6^r + 63“-/ 

-0"*031 


No other change or additional coefficient has been greater than 0"’010. 

The Addendum containing the results obtained by adding together terms 
of the same argument in Section vi. is also new. 

During the summer of last year Professor Newcomb’s new work* on 
this subject was published. His methods differ so completely from those 
given here that no comparison is made easily except in a few of the final 
results where the indirect action is insensible or is separated from the direct 
action. For the large inequality due to Venus he obtains a coefficient of 
14''’83 while mine is 14'’’55 ; a portion of the difference is probably due to 
certain terms of the second order relative to the ratios of the masses of Venus 
and of the earth to the sun which Professor Newcomb has included. On the 
other hand, he states that the possible errors arising in his method may be of 
the order of this difference, while such errors are excluded from my result. 
His results and mine for the annual mean motions of the perigee and node 
agree within 0"‘01, which is the limit of accuracy to which I have obtained 
these quantities. 

* Investigation of Inequalities in the Motion of the Moon produced by the Action of the 
Planets.’’ Carnegie Institution, Publication 72, Washington, D.C., June, 1907. 


New Haven, Conn., U.S.A. 
1908 March 27. 


E. W. B. 



GENERAL NOTATION. 


The axes of y are taken in the ecliptic of 1850‘0 and the centre of mass 
of the earth and the moon is supposed to describe a fixed ellipse around the sun 
in this plane. As it is more convenient to use the motion of the centre of 
the earth than of this centre of mass, a slight well-known change, noted 
below, is necessary in the disturbing function. 

The axis of x is parallel to the line joining the earth and the sun. 

In the scheme of notation which follows, two sets of constants are given 
for the mean distance, eccentricity and sine of half the inclination of the 
moon’s orbit. The first set is that which I have used in the expressions for 
the rectangular coordinates of the moon ; the second set is that of Delaunay 
in the final form which he gives to the expressions for the longitude, latitude 
and parallax. The longitudes of a planet , of its perigee and of its node 
are as usual reckoned along the ecliptic to its node and then along the 
orbit. 



Moon 

Earth 

Planet 

True long. 

V 

F' 

V 

Mean long. 

w^-l g + h 

T 

P 

Mean anom. 

w^ — w^-l 

V - T- w' 

V -P — w' 

Mean long, of node 

w^ = h 

0 

h" 




dP 

Mean motion 

n 

n' 

dt 

Mean distance 

a, a 

a' 

a" 

Eccentricity 

e, e 

e! 

e!' 

Sine half inclin. 

h, r 

0 


Cool’S., origin earth 

£C, y, Zf r 


rt, 1 A 

„ » sun 


aj', y\ 

y'\ »•' 


n — mean motion of the moon == 


dw^ 
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GENERAL NOTATION 


Cl, Ca, Cg are the canonical constants complementary to after the 

problem of the moon’s motion as disturbed by the sun, supposed to move in 
a fixed elliptic orbit, has been solved. 

R = the disturbing function of this problem, arising from the direct 
attraction of a planet. 

The symbols for the mean longitudes of the planets are: Mercury, Q; 
Venus*, V; Mars, M; Jupiter, J; Saturn, S. 

* No confusion will be caused by the use of the same symbol for the true longitude of the 
moon and the mean longitude of Venus. The notation of Badau has been adopted with a few 
changes. 



SECTION 1. 


THE EQUATIONS OF VARIATIONS. 


Let Wi,W2y represent the mean longitudes of the moon, of its perigee 
and of its node, and suppose the problem of the moon as disturbed by the 
sun, has been solved. Then it is well known that if a disturbing function R 
be added to the force function of the moon’s motion, the change in the latter 
due to R can be obtained by solving the equations 


/n dc{ _dR dwi 
^ dt dwi^ dt 


dR 

dci 


+ hi 


(^* = 1 , 2 , 3 ), 


where 6i, 63 are the mean motions of the moon, of its perigee and of its 

node, and Ci, C2, C3 are the canonical constants corresponding to WiyW^y 
the Ci are functions of the arbitrary constants /i, e, k of the moon’s motion, 
and they also contain the constants n', e', depending on the sun’s motion. 
The substitution of the new values of c^, Wi, thus found, in the expressions 
for the coordinates will give the disturbed position of the moon at any 
time. 


The constant part of R only gives constant additions to the hi, i.e., to 
the mean motions*: this part will be neglected, since it has no effect jon 
the new terms to be found. Hence hi = n. 


Change to the semi-canonical system n, 
changed. Putting 

dci 


dn 


= - a% 


C3, retaining the Wi un- 


and remembering that 


(la) 


dci dhf^ d/Ci dh^ db^ dh^ 

dc 2 dti dc^ difh dc^ dc^ 


* I have found these changes in an earlier paper : Tram. Amer. Math. Soc., Vol. v. pp. 279 — 
284. A fresh computation just made gives 2"*69, - 1''*42, for the mean motions of the perigee 
and node respectively. The former is 0''*03 more than the value given in the paper. 

1—2 
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[sect. 


we obtain equations (1) in the semi-canonical form 


dn 

1 ( 

^dR 

db^ 

dc^ 

db^ 

dcjN 

dwi 

1 dR 



dt 

"a’/SV 

dvji 

dn 

Tt ■ 

dn 

' d^j’ 

’ dt 

dn'^ 



dd dR 
dt dw^ 






dw2 

dt 

= + + 

(5 -‘■I 

j d^a 

' d7^’ 

dc^ 

JR 






dws 

dR , 

fdwi , \ 

dbs 

M'' 

~dw. 





> 

dt 



1 — ? • 
dn ^ 


in these equations b^, are supposed to be expressed in terms of n, Cj, 
and Ji in terms of n, Ca, Cg, Wa, w^. 

Consider any periodic term of R : 

R = n^fk^A cos (qt + q) = cos {i^w^ + Wj + (q't + 

where a is the linear constant of Hill’s variational orbit and of my lunar 
theory and -4 is a numerical coefficient (that is, its dimensions with respect 
to time, space and mass are zero) ; q't -f q" is a combination of the solar and 
planetary arguments. Then since qt-^q' is independent of the d and A of 
the Wi, the first three of equations (2) become 


dn 

dt 


'0 


a*^ 


A 


dq 

dn 


sin (qt + q), 


dd 

dt 


= - 1^71'^ A sin (qt + q\ 


dd 

dt 


hn^sJ^A sin (qt + q'). 


It will now be supposed that R contains a small factor whose square may 
be neglected. The coefficients in the right-hand members of the last set will 

then be constants and we can integrate. Put and let Bn, Bd, Bd 

denote the increments of n, d, Cg, due to R. Then 


m a^ dq n' . . 




Again, if we put 




Aj = - a^n 
48 = -a*a 


dA 

dCi ' 

dc/ 
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the other three of equations (2) become 

dwT, ^ . A . 7 


(8a) 


n “ ^ f A 1 


dw^ n“ a' 
dt 

diUs a“ 


\dt -n + y + + 

Let hhi, Bhiy Sb^, Swi , denote increments due to R. Then 

8 bi = 871 = — ^ \ ^ ^ cos (qt + ^0» 

^ dn q 'i ' 

5.1 ^6-2 5. db >2 js fj 

060 = ~i' 0?i' "t* “7 — OC2*t* J ■ 0C3 

d)i dc^ dc^ 

, n' / 1 do . .db^ . .db^ .\ 

= " + “ + “ + 

by equations (la). 

do „ do 

Putting </2 = — > ^3 = ~ ^ > 

and using the last of equations (1 a), we find 


662 = — li' 


<1 

a- 

db.i 

dn 

■ s--: 

1 cos (qt + q ) ; 

Bbs = — n' 

a- 

d-^ 

9 

is: 

db.A 

dll 


cos (qt 4- q'). 


The undisturbed values of 61, 62, 63 are of course equal to those of 

d^ d^ d^ respectively. To obtain h^Vi, bw^, (the increments of 
dt ^ dt * dt ' ^ ^ 

Wij Wz, W3), it is necessary to replace Wi by Wi-\-Bwi, and hi by bi+Bbi in 
(3a), to substitute for Bbi the values just obtained, and then to integrate. 
These operations give 


( 4 ) 




The equations ( 3 ), ( 4 ) constitute the solution of the problem. 
The form in which a periodic term of R arises (see Section II.) is 

^ = 1 ??L n'*a*il cos (qt + q'). 

4 m 
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[sect. 


It is therefore necesvsary to multiply the right-hand members of (3), (4) 
by 

Next, let 

s = no. of seconds in the daily m^an motion of the sun = 3548"*19, 


>} 

Then 

Also, put for brevity 


argument qt + q. 


(5) 1 


n ^8 
q~~s' 


/= 1 ^ ^ 206265, 

4 m /3 

l^'aW 206265. 
4 m yS 


The coefficients of the right-hand members being thus expressed in seconds 
of arc, equations (3), (4) become 

Bn ^,dq A / . , /V 
- = *,/9/ -cos(.y< + 9), 


(6) 


W' j cos + 9 ), 


(/3A+ + 5 . J)|sin(^« + ,/), 

where, to recall certain definitions, 

/S = - ^ . 9 = + I'sfts + 9", 


, rfo 


g, = - ««• 


dco 


Numerical form of the equations of variations. 

It remains to be seen how these quantities may be put into numerical 
a* . . 

form. In /, f the factor — is immediately obtained froin Hill’s results* for 

w!* 

the variational orbit ; m is a well known quantity ; — j is known as soon as 


Amer. Jour. Math., Vol. x. p. 249. 
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the particular planet is chosen ; thus /, /' remain the same for a given 
planet and ///' for all planets. The coefficients A, Ai, A^, A^ will be 
found later on, while s is known as soon as the particular term of R has been 
chosen. 


There remain for calculation 

dci dhi dbs dha dh^ ^ dh^ dhs 
^ ^ dn' dn' dn* dc^ ’ dcs ” dc^ * dc^ * 

which, depending only on the orbit of the moon as attracted by the sun and 
earth, are the same for every perturbation of the moon’s motion, and there- 
fore apply not only to the present investigation but also to all investigations 
where a disturbing function R is added to the moon’s force function. 


Some idea of the degree of accuracy required is desirable. The largest 
known inequality is that with the argument which has a 

coefficient of about 15". For this = 1, i.^ = — 1, ig = 0. The principal part 


is given by 


.A dq 

d 8^ dn^ d ' 


dh.\_ 
dn) ~ ^ \ 


-(1+ 01486). 


There is no other coefficient which is so great as 2". Since the degree of 
accuracy aimed at is 0"’01, it will be sufficient to use four place logarithms 
and four significant figures for the functions (7) so that the final results 
will be accurate to at least three significant figures. But certain of the 
functions are only needed to one or two significant figures, as will appear 
immediately. 

The functions Ci, C2, C3 are the same as Delaunay s L, G — L, H — G after 
the final transformations and the changes to his final system of arbitraries, 
n, e, 7 have been made. As my results will be used for the calculation of 
the moon functions, it will be more convenient to transfer A^, A^ to my 
constants e, k. 


Let ^ denote the derivative of a function Q with respect to n when it is 
dn 

expressed in terms of n, Ca, C3, and when it is expressed in terms of 

n, 7*. Then the following equations serve for the transformation of the 
derivatives of Q from one set to the other*: 



The functions considered here involve y only in the even powers. 



8 


THE EQUATIONS OP VARIATIONS [SECT. 


the same equations will serve for the transformation from the set n, Cg, Cg to 
the set 71 , e®, k* if we replace e®, y* by e*, respectively. 

It is first to be noticed that 


Cj = naV ]— J 4- power series in tti, 


Ca = TittV {— 2 -f 


Ci — Tia^ { 1 4 »* « „ » } > 

{ J-l- „ „ „ „ }, 

&3 = nm' (-|4 „ „ „ „ ); 

the power series in each case vanishing with m. Hence when for Q in 

equations (8) is put Cx,h^ or 63, the principal part of the first term in each 

equation will have a portion independent of 7®, while the other terms will 

have one of these quantities as a factor. Since 3^, 7 — A approximately, 

two significant figures will be fully sufficient for the values of . It is 

an an 

to be noted that the second and third of equations (8) do not depend on 
derivatives with respect to n, and therefore that we may change to the variables 
e, 7 or to e, k, from Cg, C3 without reference to the first equation. 

From Newcomb’s transformation* of Delaunays values for L, (r, H, 
we have 


Terms in 

Tui^e^ 

A. 

aV“ \dn ) 

Terms in 

r.s 

na'iy^ 

a'^'^ \dn J 


-•50015 

+ -1667 

7)1^ 

-1*99704 

+ •6657 


+01686 

- -0393 

7rr 

- -00322 

+ -0075 

7n^ 

+ -00644 

-•0215 

7n^ 

+ -00049 

-•0016 

7)1* 

+ 00170 

- 0074 

7)1* 

+ -00005 

- -0002 

rem. 

+ -00032 

- -0020 

rem. 

- -00022 

+ -0010 

Sum 

-•47483 

+ -0965 

j Sum 

- 1-99994 

+ •6724 


The remainders in the first and third columns are obtained from the 
values of Cg, Cg (calculated by processes independent of the nature of con- 
vergence along powers of m) given by myself f ; those in the second and 
fourth columns are estimated from them. In any case the results for 

(^) ’ (dn) correct within two per cent, and this is all the accuracy 
necessary for our purpose. 

In a similar manner the values of can be obtained with 

* “Action of the Planets on the Moon,” Amer. Eph. Papers^ Vol. v. Pt 3, pp. 201, 202. 
t “Theory of the Motion of the Moon,” Tram. R. A. S., Vol. lvii. pp. 64, 66. 
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more than needful accuracy from the series in powers of w given by 
Delaunay*, Hillf, and Adams t with my numerical values §. The deri- 
vatives with respect to 7 * are obtained immediately from the last named 
reference. The derivatives with respect to e^ k“ can be immediately derived 

by inserting the values ^ which I have also giveii§: the change in the 
derivatives with respect to n is insensible. 

I find the following results to four places in the logarithms, to each of 
which 10 has been added : 

= + - [9 0775] na^ - [7'5002] 

(^) ^ a^ + [6-7577] na\ = - [10-3039]«o=, 

(a)-!*-'™*!’ S— »-[»’“«’■ 

(g)- + [7-578«), + 

Whence, accurately to four figures, 

-[81720], +[#817.5], «.jJl. + [83»0J, 

1;- + [7.5738], .•'^. + [#■3900], [7.5C98), 

We have also 

^^=+[999921], m = + [8-87391], /3 = + [9-51801], 

' /= [22-29358] /' = [17-6174.8] = [5-32390], 

/3 being found accurately in Section ill. From these quantities the terms 
depending on - in (4) can be found. 

Further, by substituting A*, successively for Q in the second and third 
of equations ( 8 ), and putting e for e, and k for 7 , we find 

J.-a.»"-[10.7.»l"-(87782]". 

* Comptes RenduSf Vol. lxxiv. pp. 19 et $qq, 

t Acta Math., Vol. vin. pp. 1—36. Annah of Math., Vol. ix. pp. 31—41. 

X Monthly Notices, R. A. S., Vol. xxxviii. pp. 43—49. 

§ hoc. cit., Vol. LXiv. p. 632. 
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Finally, making these various substitutions in (3), (4), the equations of 
variations become 

' Sn A 

^ = (— 4- -0148612 - *0037 444)/' — cos (qt -f 9'), 

= + [9-51801] hf ^ cos {qt + q'), 

= + [9-51801] 4/' j cos {qt + q) ; 

Bw, = |(- i, + -014804 - 0037444) /- + /' ~ 4 sin {qt + o'), 

(9) J I s s ) 

Sw, = |(4- -014864 - -0070664 - -0081484) 

+ (- [8-1720] A, + [11-0999] ^ - [7 9422] -^'J sin {qi + q'), 

Sw, = |(- -0037444 - -0081484+ -0012104)/^ 

+ (+ [7-5733] A, + [10-2620]^^- - [82962] ^)-~| sin (qt + q), 
where A,^ — ^ (As?), 

a^ (171 ^ 

A being expressed in terms of Ca, the figures enclosed in square brackets 
being, as elsewhere, logarithms with 10 added. 

These equations replace the equations of variations for Bl, Bg, Bh, 3a, Be, By 
first given by G. W. Hill* and calculated also by Radauf. To reduce the 
above system to theirs it is necessary to find 3c2, Bc^ in terms of Bn, Be, By. 
But it is simple to compare the first terms of Bwi, on which the principal, part 
of each long-period inequality depends. Radau finds 

Bwi = {(— 3‘0676i 4- '05601i' — *01124i'')p -f ...) Psin {qt + /), 

where 4 = i + i' + i", ^ = + 4 = 4', = 

This reduction gives 

= |(- -98944'+ -014764 - -0037054) / ^ + . . . | sin (ji + q'), 

which is about less than my value. But I have been unable to find the 
coefficient — 3*0576i in Bwi, which Radau gives: from his data I make this 
coefficient —3*0791]; which, reduced to my form, becomes — *9964. The 

* American Eph. Papers, Vol. ui. p. 390. 
t Loc. cit. , pp. 86, 86. 

This apparent error seemg to be due to some confusion in the substitution of the numerical 
values for n before and after the final Delaunay transformation. 
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difference, *0036, is about the error to be expected owing to the slow con- 
vergence of the series from which Radau obtains his coefficient. It may be 
added that the value of on which my coefficient of ii alone depends in this 
connection, can be determined accurately to at least seven significant figures, 
and the value for it quoted aboVe has been again verified and tested. 


The Variation of the Moon's True Longitude. 


The coefficients in longitude are obtained by substituting the values thus 
found in 

(10, 8 r.o., + (g-i)s., 

dV i ,dV . dV . dF . dV ^ 


the first term of which contains the primary inequalities, and the remaining 
terms the accompanying secondary inequalities. 


The largest coefficient in V — that of sin^ — produces a maximum 

dV 

coefficient in SK through . Sco, less than 1"; the coefficient of sin ^ in 


V is 22640". It has therefore been necessary, in order to include all 
coefficients greater than 0" 005, to find those terms in V having the factor 
which have coefficients greater than 100", or, if jo = 0, to 400". 


From my results I obtain 


[9-0405] sin / -f [8-3470] sin (2D - 1) 

+ [0-9688] sin (2D + 0 + [6-2703] sin (4D - 1) 

+ [8-0603] sin 2D + [7*57 15] sin 21 

+ [7 0112] sin (2D - 21) + [6*1738] sin (4D - 20 
~ [7*3001] sin 2F + [6*4273] sin (2D - 2F) - [6*7818] sin D 

- [6*7262] sin (I + 1') + [7*0000] sin (2D - ^ - 1') + [6*8555] sin {I - 1') 

- [6*3398] sin (I + 2F) + [6*2826] sin {I - 2F) 

+ [6-2434] sill 3^ + [5 -8060] sin (2D - SI) 

~ [7*5110] sin V 4- [6*9040] sin (2D - V\ 

the coefficients being expressed in radians and the notation for the arguments 
being that of Delaunay. From these results the derivatives with respect to 
Wi are immediately obtainable with 

D = Wj — earth’s mean longitude, 

F = w;i — Wi. 



12 


THE EQUATIONS OF VARIATIONS 


[sect. 


For ^ the results of Delaunay might be used, but it will be shown 

dV 

immediately that hn contributes nothing sensible to SF — a result 

probably of the use of the system n, Ca, Cs instead of n, e, 7. The derivative 
has however been found for the largest terms in order to show this fact. 
The derivatives with respect to Ca, C3 are obtained from my results with the 
help of equations (8). They give 
dV 

n + [8-0465] sin I - [8-4669] sin (2D -1)- [7*2838] sin (2D + 1) 

- [8-3965] sin 2D + [6*8784] sin 2/ - [7*1761] sin (2D - 2^ 

-[6*8271] sin 2F, 
dV 

na^ - [11-5828] sin I - [10*8895] sin (21) -1)- [9*5059] sin (2D + 1) 

- [10 0018] sin 2D - [10*4147] sin 2/ - [9*8546] sin (2D - 2Z), 

= [9*6959] sin 2F. 

The latitude and parallax are treated in a similar manner. 

Abbreviation of the Formulae for 3F 

In the actual applications to the calculation of the direct inequalities 
much abbreviation of these formulae is possible. For example, in the large 
class which has its arguments independent of the lunar angles, = i^ = i^ = 0, 
and therefore hi = 8ca = = 0. The maximum values of the remaining 

coefficients are 

/'M- 

s ’ s de 5 ^ s rfk id ’ 

. dy_. 1. ^ 1. • A 

dwi ’ 9 ’ dw,2 ’ 9 ’ dws ’ 250 ' 

Hence, to the accuracy desired, we can put 

sin (^^ + q\ Bw^ = [ 1 1 *0999] ^ sin (qt + q), 

= [10'2620] + q'). 

When t’a, 1*3 are not all zero, it is still possible to limit the formulae 
very materially after Bwi has been found, owing to the limitations on the 
magnitude of the coefficients in hv^ and F. 

For this purpose we note that, for the greatest coefficients, 

\dV\ ^ 1 \dV\ ^ \dV\ 1 


1 

8 dk id ’ 

• "Z Z 

dw,’ 250' 


dV 1 


AdV 

dn ^30’ 

, .A 

wa* j— < 40, 

acj 

Idc, 


-1 


dw^ I 250 ' 
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Put -f -014861; - *003744t;, 

^ ~ S»” ’ 

= [5-32390]. 

^ ^ = [5-3239] Cs cos {gt + g'), 

Then (11) ^ J^4 = [4-8-«9 ]^*^cos(5< + ,/), 

^ [-*'8419] i,cos (qt + }'). 

Also, if Xj = -01 +86t. - OOTOeO?; - -008 1 481, , 

- 0037441, - OOSHSf,, + 0012104, 
dA 


^ “ dk 


, dA . , 


the equations for hwt become 

Sw, = ^ • ^') (9* + 9'). 

8w, = J- [3-4959] J + [7-3842] 

- [4-<>160] 6' sin {qt + q'), 

8w, = j+ [2-8972] ^ ^ 4- [.5-9349] j. 

- [4-5805] *- j, + M 0 sin (qt + q'). 

Ai) 

Now C<lo", ■^^' = 0^®-.^'<0"-1. j <1, A. 5 2. 
s j A Ai 

Hence < 5000", Cs < lOOOO". 

Ai 

Therefore if the limits of the coefficients to be considered be 0"'003, 
dV 

^ Sn and the parts contributed by the first and third terms of Sw^ and 

hwt are insensible ; as a matter of fact no one of these parts is so great as 
O^'OOl. Hence we obtain the final form of the equations of variations given 
on the next page. It is necessary to remember that, although the derivatives 
with respect to c,, c, have been transformed into derivatives with respect to 
e, k, the derivative with respect to n is to be taken on the assumption that 
the coefficients are expressed in terms of n, Cj, c,. 
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[sect. 


Final form of the abbreviated Equations of Variations. 

+‘^') ““ (9i + q) 

= c(l+ [5-3239] ~ . sin {qt + 5 '), 

Sws = |[7-3840] ^ sin (qt + q’). 


^,= |[4-8419]^_4C'Jco3(9< + g'), 
hw, = |[6-9349] ^ Csin {qt + q), 


where tt = \ n^a.’A cos (at + o'), 

4 m Ml 

X, = _ + -0148013 - -0037444, 

\3 = + -014864 - -00706G4 - -0081484, 
X3 = _ -0037444 - -0081484 + 0012104, 

4> 4> 4 'ii’o the coefficients of w,, w^, w, in qt + q', 


“ ni A., 


(lA . 


j,A, k =j,A. 


One or two particular cases of frequent occurrence may be mentioned. 

fA 

In all cases 0 — . For multiples of the arguments given, divide the 

second term in Bwi by the multiple. The argument is the moon portion of 
qt + q. The omitted terms are either zero or negligible. 

A rg. 0. Here ii = i^ = i^ = 0, 

f'A 

8wi—' ~ ^ HUi (qt -i- q\ 


1 A 

[110999] j ~ 8ivi, 

Sw, = [10-2021] Sw,. 

Arg. J = w, — w,. 

Sw, = |l - [5-3175] s ^'1 Csin (qt + q), 

Sw, = {- [7-3776] s - -02161) G sin {qt + q’), 

^ = [4-8355] Cs cos {qt + q'), 

StVs = — [6 9285] sin (qt + q') (for the latitude only). 
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Arg, 2D - ^ + W 2 2^. 

St(;i= |l - 

Sw 2 = - [7*3905] sG sin (qt -f q), 
- [4*8484] sO cos (qt + q'). 


[5-3304] s 4^1 C sin (g« + g'). 

A 


Arg. 2^ — 2 D = — 2 w 2 + 2T. 

Swi= |l - 

Sw 2 = {- [9*21 19] s- *4753} Osin (qt + q), 

^ = [6-6698] sC'cos(g«4-9')- 

There is quite an extensive class of terms containing this argument in 
which Swj is insensible. For these we can put 

Swi = [5*7642] *- . ^2 sin (qt + q'), 


[6-8505] s sin (fy« + 5 '). 


^^J=-[3-2221]-{.^cos(9t + g'). 

The great majority of terms in R to be considered contain powers of the 
lunar eccentricity as a factor and the principal terms in the moon’s true 
longitude have the same property. In nearly all such cases it is permissible 
to neglect higher powers of the lunar eccentricity and inclination. When 

these two conditions are satisfied it is not necessary to find Sc.^, ^ . l?or 

then jo = Uq 1 and the ratio of the first term of the coefficient in Sw^ to 

. r I n. • . dV dV rpu’ 

that of Scj is the same as the ratio of the coefficients in c . Inis 

arises from the fact that Sw 2 depends mainly on e ^ , while Sc^ depends on 

dR 


dwo ' 


If then 


= Q sin (ij Wfl + F = O' sin (i^'w^ + ^|r') , 

where yjr, are independent of «/a» and t’a, tg' have the same sign, 

^ 8w2 + ^ s= QQ'h' sin {(4 - iV) Wg + ^ - 'vfr'}. 
dw^ ofCji 

It is to be noted that Q is here the first term in the coefficient of I 
An exactly similar theorem holds with reference to c,, and the terms 
which contain powers of k as a fector. 
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THE TRANSFORMATION OF THE DISTURBING FUNCTION. 


It is well known that the disturbing function for the direct effect is R, 
where 


R _ 1 + y r] -F 


In order to take into account the motion of the earth round its centre of 
mass, the terms containing are as usual to be multiplied by the ratio of 
the difference and sum of the masses of the earth and moon. 


Put 

then, by Taylors theorem, 


d d d 


_a_ 

dQ’ 


(13) 


1 =ri_i + iy9v_i/AY+ 11 

[(I _ a;)^ + _ rjf + (f- L 3Q i2 VQ! |3 V0Q/ J A 

1 a;f + yTj + rf r 1 / 9 y 11 

“A+~ A» ■■■Ja' 


Since R will enter only through its derivatives with respect to the lunar 
elements which arc not present in A, we obtain 



The separation of the terms of R into sums of products, one factor 
involving the lunar coordinates and the other the planet’s coordinates, 
suggested by Hill* and adopted by Radauf, is implicitly used here, and the 
division is the same in principle, but it will take a different form with the 
use of complex coordinates, and the method of sepaiation is made on a 
general plan which can be applied to terms of any order in R. 


* Jmer. £ph. Pap^rg, Vol, in. 
t de VObi. de Paris (M4m.), Vol. xxi., 1892. 
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aj-l-yV— l=u, f + l = Ui, 

x — y V- 1 = s, ^ — V 1 = Si ; 

then A’* = UjSi + r’' = us + 

0 0.0 0 
0Q ^0Ui ^ 0Si ^0?‘ 

And since 1/A is a solution of Laplace’s equation with respect to f, f, 

a?” A “ laf aijv A Sihas, a • 

Hence, omitting 1/A for the sake of brevity, 

1 / a \» 1 / 3 V , 1 / 3 \V S'" 

2 (aQ) “ 2 (“auj 2 V as,) an.'as. 


l,/3\\ 3/ 3, a\ 

2 ^ (s?) ^ 3^ ('' 3u, ^ ® BsJ 


= real part of (»- - j + 2^ 

' = real part of - (»- - 5.^) n g^- - J (g^ ^ 




8 /0^011,08/ 
/ 0 \ , 


( 0 \ 

0(^/ replace 

\ 27 n / pi-2 \m 

3?) 

I now give a method for replacing the derivatives with respect to the 
coordinates, by derivatives with respect to a', 7 ’ (the earths moan longitude) 
h",y". In this way it is possible to make use of Leverrier’s expansion* o 
1 /A in terms of the elliptic elements of two planets to the seventh ordei 
inclusive with respect to their eccentricities and their mutual inclination. 

Expressing the coordinate.s in polar coordinates, we have 

/ q* (1 — t'" cos ( V" — + y'‘^T'' cos + P — 2/t ), 

( 15 ) - ^ = (1 - 7"") sin ( V" - P') - y'^r' sin ( P" 4 - P' - 2 ^"). 

X = 27" Vr=VV'' sin ( P" - r). 

Whence A’ = r"^ - r * - 2 /f , 

Ui = — r + (1 - y'^) exp. (P" — P') V— 1 

q- yV exp. (- P" ~ P' + 2^) V-T, 

the expression for Si being obtained by changing the sign of V— 1. 

* See below. 

BR. ^ 
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From these results we obtain 

d _d 

dr' 0u, 0Si ’ 


V-l A = 


(16) =r + v~l ^177= u,r Sjx 2r 

'0v 0r 0F 0Ui 0Si 0Si 

0 , 0 / z 0 9 0 ^,0 

9 0 

the notation , g— being introduced for a few pages for brevity. 
Let F denote a function of ^ and of the product UpSi only. Then 


(16a) 

(17) 

Hence 

(I.-) 


dF ^ ,dF dF ^ ,dF 

0V 0Si ’ 0W ^ 0Ui ’ 

/ 0 1 _ / 0 L — u " F 

[duj [dj A“^'^ 


0 = 1 _ 
0V0W ‘ A 


(r' 9,y + (AYli = l/ 2.'^ 

\ dr'J \0F7 J a 0v \ 0 


d_ 

9 uj 


= 4 r 


d'^ 


d_ 

0V 


1 

A 


\9v 

= 


0Ui0Si * A ’ 


02 1 


Similarly, 


0s, 2 A * 


1 1- i.(- 2 rF-^-" 

0v 0w A 0Ui0Si2 A ’ 


and the general law is evident. The derivatives with respect to u,, s, are 
thus e.xpressed in terms of derivatives with respect to V\ when for 
0 0 

^ have been substituted their values (16). 

0V 0w ^ 

Next, 

A^ = r'2 + — 2r'r' (1 — y '^) cos ( V" — V') — 2rr"y"^ cos ( F" + F' — 2h") 

= r'2 4- r"* — 2rr' cos ( F" — F') + ^rr"y"^ sin ( V' ~ h") sin (F' — h"). 

1 _ 2 (1 - 7 "*) rr" sin ( F" - h") cos ( F' ~ h") 

A A> 

0 1 2(1- 7 "*) rV" sin ( V" - h") sin ( V' - h") 

07">A A» 


Hence . gj,) 


(1-7"*) 
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Therefore 


0 1 d j — :: 


' ' ’ ^ > dy''^} A 

2 (1 - 7"0 r'r" sin ( 7" - h") exp. - ( V' - h") 




exp. -V- 1(7' -A"). 


7" A’ '■ 

This gives 

0 !■_?’_ 7 '' ^'^P* V— i ( V'' — /i") 

0f * A ~ A^ Vl — 7^2 r 


1_ ^ A 

1 17' ■!“ 6*. "2 


+ V-l(l-7-)g|.,}^i. 


10 7'^27"^0r ■ " "07'^J A ■ 

(As f/A® is real, we incidentally have the relation 

[.in ( V - r) . - C ( r - r) (1 - y.) ,[,„] I . o, 

and the sign of the V— 1 may be changed.) 

Since f/A'^ involves Ui, s^ in their product only, the formulae (IGa) are 
available ; and as 

0 exp. \/- 1 ( F' - h") o_ exp. +y- j ( ^ - 2^ F 


0 exp. \/-l(F'— ///') ^_ex p. 4- V— 1 (F' - h") ^ ^ 
0w r r' 0VV 

we obtain, by substitution of the various results, 

( 1 ») 


P- 5., (jf- 


9 . 1 9 

)wl97''^27"“9r 


+ V-'l (1-7"'') 


97"“J a 




+ [" ("*-§"“) ( 1 ^ - 9 -^ 

9w] [gF + 2 y^> 3 F' "P ^ e - 7 ' ’) g”^i] 
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The derivatives with respect to v, w are easily expressed in terms of deriva- 
tives with respect to the elements by means of the well known formulae 

'1- ' L A.- i 

dV'~'dT’ 

where T is the mean longitude of the earth and where it is supposed that the 
angles referring to the earth are expressed in terms of T, I'. Also, Leverrier’s 

development depends on a = — , where is the mean distance of the inner 

(I2 

planet, that of the outer. Hence if we put 

0 


da 


1 


we have, since is a homogeneous function of degree — 1 with respect to 


Ct f d 





0 

di — a , 

II 

dr ~ 

^0a 




0 

e 

il 

a .2 = a'\ 

dr ■" 

“0 a 


— 1 = — i) — 1 for inner planets, 
= D for outer planets. 


The 'planet portions of the disturbing function are therefore expressed in 
terms of derivatives of 1/A ivith respect to a, T, //', 7 "^^, — the quantities present 
explicitly in Leverriers development. 

Each term in the development of R consists of two main parts together with 
a constant factor. The first part is a function of the moon’s coordinates alone 

^witli which will be combined the factors ; the second 

part consists of the derivatives of 1/A which is a function of the earth’s 
and of the planet’s elements only. Let 0 denote an angle present in 
the first part, an angle present in the second part, that is, in 1/A, so that 
^ ± </) is the argument of a term in R which is under consideration. Put 
j for V~l, and* 

7^ * — a® — ~ 

%■% = i (M + M,) ei> M,) e-^ = M, cos 0 + jM, sin 0, 

% ^ ^ 

These functions are calculated once for all and serve for all the planets. 


* The exponent always attached here to the exponential e prevents any confusion with the 
same symbol used to represent the lunar eccentricity as defined by Delaunay, 
t Owing to the absence of h" from z, u, /, the part is not present here. 
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For the second part put* 

i = P cos (j> = P cos {iT 4- ^iih" + 0'), 
where (f> is independent of T, h". Also for inner planets 

• !:=[(»•' = + 

2P, cos ,!> - 2 P 3 j sin - 2 ) ^ 

= (Z) + 1 + i A) (d + 3 +;• v) P cos i>, 

Pt cos <f> + jP, sin 4>=-(j)+l+j 


so that 


P, = [(P + l)>i-t‘‘] P, 

P, = i[(D + l){D + ii) + P]P = iP: + iI)+2)P + {P-l)P. 


P, = (D + 2}iP, 

P, + P. = (P + 1 - i) J ^ + + (1 - 7 "“) A 

P, - P. = (P + 1 + i) {-*•- 4, + (1 - 7"“)g-|4 P- 

The first line of P, which consists of tViose terms in (18) which have the 

factor , becomes with these substitutions, after inserting the unit factor*|' 

n^OL^ 
m ’ 

(19) ^ \m,P, + M^P^ T M>P, 

- ± ± 

either all the upper signs or all the lower signs being taken according as it 
is convenient to use the sum or difference of the angles 6, It will be 
noticed that the Pp contain the divisor a which will cancel the a outside 
the parentheses; otherwise they contain a\ a" only in the form aja"=^a. 


* The letter i has also been used elsewhere in a different connection ; but no confusion need 
occur. 

t The denominator is properly the sum of the masses of the earth, moon and sun ; but the 
difference is insensible. 
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For exterior planets the same formulae hold if (1) we substitute — i) — 1 
for D, (2) replace the factor a by aa'\ since Leverrier’s formulae are then to 

be expressed in the form func. (a). 

8» 

For the second part of R, which involves the factor ^ and which will be 
denoted by R^, put for the moon portions, 

u (r» - 5 ^») ^ ^ ^ (Jtf, _ M,) e-^i = M, cos d sin 6 , 


r a 

1 ^ 

8 * a^ 

a 


= ^ {Ms + Ms) i (^8 — -^9) = Mg cos d + jM^ sin 6 , 


iz 

'3 • “ — ^ ““ ~ = -^10 (cos 6 +j sin 0 ), 


r'^ ' a‘* 

and for the planet portions, in the case of the inner planets, 


= = ■^12 (cos 6 + j sin 6 ), 


P„ cos <l> -jP, sin <^ = - 2) ^ . P cos </) 


= _ (z) + 3 + j A') [(Z) +\y- i»] p cos 4 ,, 


P« = {dt - (i " 0w 

= -^P + 5+j^,) (P+3+i 


. d 
dV" 


P„ cos 4 + ji’„ sill </) = j - 2 ^ 


I 


a 

0W 


^P+l +jgy7) PCOS^, 


0 r 27 "= 0 A" ■'■•^ 0 7 ') 0 y>a[ -P 4 > 
. 0 




. f 07= + (0^ + 27 "“ 0 f)} ^ 

P„ COS 4 > + jP., sin 4 >=j - 2 ) , 


0W 


■^ + 4.,+i(l-7"n0^jPcos^ 


= (P + 3 + j 0!^) (-D + 1 + j gir-) 


|- (1 - 7"’) 0^--, +3 (0 + 2y'. 0/,")| -P 4 >. 
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SO that 

■P. + P 7 = - + 1)’ - (i) + 3 ± t) P, 

P, = -(P'“ + 6P + 5 + 3P)(i)+3)P, 

P, = - iW‘ + 18 P + 23 + tO iP, 

A. ± P„ = - {(P + 3) (P + 1) - T 2*} {(1 - 7"^) + {i + P, 

Pn + P.a= -KP +3)(P + l)+*n 2i(P +2)1 {(1 - 7"=) (»■+ ^ 7)1 P- 

Then 

(20) = ~ a' | JI/,P. + M,P, + M,P, + ilf.P. 

ib m a 

- {Mr,iP,, ± Pn)+M,,(P,, ±P„)) j COS {0 ± <^). 

The changes for outer planets are the same as in tlie case of Ri. 

The moon coordinate u is referred to the true place of the sun. As my 
results will be used to calculate the functions of u, s, z, and as these results 
are referred to the riiean place of the sun, it is necessary to replace u by 
where Uo is the same as the u of my Lunar Theory. 

The formulae for Mp and Pp are collected and slightly altered in form at 
the end of this Section. In appearance they are somewhat complicated, but 
they arc easy and rapid for numerical calculation. For example, the effect 
of the first factor in P4, P5, Pjo, Pn, P12, P13 on the expansion of 1/A is 
seen at a glance ; it is shown below that where Pj has not been tabulated it 
is rapidly found from a general formula; and the effect of D on P or P, 
requires nothing more than the replacing of Kp by (/> 4- \) Kp^i-\- pKp, The 
great majority of terms involve only Pi, P2, P3, while the functions 
Pg, P7, ... are required only for the very few sensible terms which have 

the factor 

a 

In the above development it is assumed that the elliptic values of f, rj, f 
are used. If the effect of their perturbations is to be included, the new part 
BR of R can be found by operating on its value (14?) with 

where 877, are the disturbed values (supposed known) of f, r}, f. 

Similarly, x, y, z represent the coordinates of the moon as disturbed by a 
sun moving in a fixed elliptic orbit. If 8t/, hz are the perturbations 
produced in a?, y, z by the value of R in (18), we can include the effect of 
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these additions to Xy y, z (perturbations of the second order relatively to the 
masses) by a further disturbing function 


but then the Wiy Ci must be newly defined, that is, they require slight 
additions — probably insensible — in the equations of variations. 


There is no difficulty in the calculation of these two new functions. 
They are obtained easily by the method given above. It will be found that 
Mly for the first term of li, involves the functions P,,, ..., P 13 instead of 
Pj, ..., while the expansion of h'R, which involves changes in the 
Mp only, is so simple that its value is obtained immediately. 


The formulae to be used are collected below. 


The Functions 



iif,= 

= coef. 

of 

ejo 

in 

a'^ r^-3z^ 
* a^ 


4-1 

M,) = 

= coef. 

of 


in 

IL g 42 j|F'-r) 
r'- 



M,= 

= coef. 

of 

e^^ 

in 

^^p}(T-h") ^ 

Ml 







7''“ 

a** ’ 


M,)- 

= coef 

of 


in 

pJjiV'-T) 

Uo (r® - Sz~) 
a=^ 

i (-^^0 ± 

M,)-. 

= coef 

of 


in 


Uo-'' 

* 3a'” 



= coef 

of 

eJe 

in 

r ’ 

zj {r^ — 
a? 



- coef 

of 

eje 

in 

(j'i 

a^ 


The Functions Pj, P^, 

^ cos ^ = SP cos iiT + 2i,A" + f ), ^ = « ^ • 

Inner Planets. 

+(Z) + 2)P + {i”-l)/', 

± A = {(1 - y"^) ^ ± (i + y‘,)} (O + 1 + i) p. 

P,±P, = -{D + Z±i)P,, 

P, + /’, = - (y)± 3t + 7) A-4(l±t) (2 + i)(i)±i+ 1)P, 

Pn ± i".. = - [( 1 - y'“) ± (^ + [P, + 2Dr + (2 + 2i) P], 

-P» ± j*., = - [( 1 - y"’) 3^, ± (*• + 4 ;.)] [ A + (2 + 2t) /jp + 2 (.• + 1 )> P], 
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Outer Planets. 


P, = \P,-DP + i^P, 

P, = -i(Z)-l)P, 

P, ± A = - (/) ± i) {(1 - r"^) ^ t (i + 4,)} P, 

P(s^P 7 = (^~2 + i) Pi, 

A±A = ('D + 3i-6)P, + 4(l±i) (2 + i) {D + i)P, 

Ao ± A. = [( 1 - y'“) ± (^ + 4 ,)] [- A + 2PP ± 2tP], 

Pn ± A, = - [(1 - y"’) ± (» + 4i)] [A + (± 2i - 2) PP + (2P + 2i) P]. 

The Disturbing Function is Ri + Ri where Ri has the value given by (19) 
on p. 21 and the value given by (20) on p. 23. 

Owing to the fact that is nearly equal to + and to + P 3 in most 
cases, the portion + M^P^ of was usually put into the form 

h (M, + M,) (A + Pa) + i (M, - M,) (A ± A) 

for computation ; the degree of accuracy of the tables in Section v. was then 
always sufficient. 



SECTION III. 

DEVELOPMENT OF THE DISTURBING FUNCTION. 

Leverriers Expansion 0/ 

Leverrier* expands in powers of the eccentricities of the two planets and 
the s(piare of the mutual inclination of their orbits. The arguments are their 
mean motions, the longitudes of their perigees measured along the fixed plane 
(ecliptic at a given date) to the nodes, then along the orbits, and the longi- 
tudes of the nodes. As one of the planets in the present investigation is 
taken to be the earth, and the fixed plane the ecliptic, supposed immovable, 
Loverrier’s arguments become in the present notation, for an inner planet, 
T = T' = r, 

X = P, ijt' = P - I". 

(Leverrier denotes the mean longitude of the earth by V \ I use this letter 
for the mean anomaly.) 

For the coeHicionts I use the .siiine notation with the exception of 7;, 
which I denote by 7 For an outer planet the accents (except in the case 
of }i \ y") interchange. 

Put Ao' = 1 + — 2a cos (P — P), 

then with Leveriier, using to denote suminatiun for integral values of i 
from x to — X , I put 

“ = ^a ' XA cos i{T — P), 

=ia':t/i''icost(r-P), 

=ia' 2 C«'cost( 7 ’-P), 

^ cos i {T — P), 

* AnnaUs de VOb*. dt ParUy Vol. i. The terms of the eighth order have been computed by 
Boquet, ih. Vol. xix. 
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and* 


or, in general, 


=io'S3'*'cost(r-P), 

cos i{T- P\ 

XV = cosi{T- P). 


If be any one of these coefficients, put 




1 

R = -- ar _ R <*» 


Then, when the functions have been calculated, wo have all the 

materials for obtaining the functions Pp. 


Leverrier, besides the notation just given, uses certain functions of 
A, B,..., in the development. The notation adopted here is the same as 
his. It is 

^ ^ -f 

QH ) ^ a (00-2) ^ 40(0 4. 0('+2)), 

(20 ^ jji, ^ I + a<'i), 

S‘‘> = If (P<‘~” + 3Z)'-" + P“+’>), 

yc-i = + D"'-"). 

Of these, it has been rarely necesjsary to use any but 6 r<o. 

As powers of the do not occur in any of the formulae which will 

be used, the brackets round the i in the index will be omitted. Fur brevity 
we shall put 

A.o' = 

and when the index i can be omitted in any equation without causing 
confusion, it will be dropped. 

The coefficients in Leverriera expansion of 1 /A are all functions of 
e\ e'\ 7 " and the Ap, Bp , .... The method outlined here for deriving all the 
planet functions from this one expansion does not necessarily give the shortest 
algebraical expressions for the coefficients, but for numerical computation, 
which is the principal end in view, these expressions have this advantage — 
that they require very little use of logarithm tables. The calculations consist 
mainly of additions, subtractions and multiplications by integers less than 
100, and the functions are read straight from Leverrier’s expansion. Moreover, 


* The reversed letters 3, are used to distinguish from other functions defined on this page; 
Leverrier does not need the two last expansions. 
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it is possible to see almost immediately when the terms in a given coefficient 
become insensible. For example, in the case of the great inequality due to 
Venus, argument Z + 16T— 18F, the principal term has the factor 7"'^ and 
argument Z + 16 r- 18 2/i", and the terms of order 7"® are just sensible. 

The following table of the coefficients actually calculated show.s immediately 
where the calculations can be stopped. The first column gives the order 
of the particular portion of the coefficient, and the remaining columns the 
parts of this order, contributed by P,, P3), iCPa-Ps), respectively. 

Argument ^+ 16 P — 18 F— W'. 


Order 


^2 2 

n>> 

y 

- 15"-89 

- l"-28* + 

y"" 

y"^ 

y"« 

+ 1*65 

+ *17 - 09 

- 11 

+ 01 

1 

6 

0 

0 


+ -27 

•00 ' - -01 


+ -12 
j - *00 

- -01 

+ -01 - -01 

Sums 

- i:V96 

- Ml + 1-08 


All of these were fully calculated with the exception of that of order 7"*, 
which a brief examination showed to be between — 0"'0()5 and — ()"’015 ; 
the terms of 01 dors e \ JGid a portion from the third term of R are 

of about the si\me size, but they were not calculated. I believe the result, 
— and the values for the other terms of this period, are accurate to 

within 0"*05 ; the adtlitional computations nece.ssary to obtain the final 
coefficient within 0"'()05 are not very long, but in view of the uncertainty 
in the mass of Venus, which is doubtful within one per cent., and of the 
length of the period of the term, the present results are fully sufficient. 

Calculation of the Functions 

The known methods used for this purpose f have been modified in order 
to avoid the want of symmetry which makes them troublesome, and the 
considerable loss of accuracy which occurs for large values of i. The point 
to be considered in this connection is not the numl)er of places of decimals 
but the number of significant figures, and the methods used here have been 

* In the original essay these portions were given to be -l' -90, +1"*16, respectively, 
owing to a numerical error. 

t See Tisserand, Mfc. Cil., Vol. i. Chap. 17 ; and Leverrier, loc. cit. 
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adopted, partly to obtain the same number of significant figures for every 
coefficient, partly to avoid numerous multiplications by incommensurable 
numbers (involving the use of logarithms) and partly to obtain easy checks 
on the results, so that there shall be no doubt about the numerical accuracy 
of the tables giving the functions It is true that functions for values 

of i not afterwards needed have been found, but I believe that the gain in 
accuracy more than counterbalances this defect. The sieve has shown the 
maximum value of i which is likely to be required in the case of each planet, 
and the tables have been formed up to this maximum value of i. The 
advantage of the present methods has been obtained mainly by the use 


( d 

a in preference 



We have the well known formula* 


( 22 ) 




2.4.6“.. 2i 




2i “h 2 


a- 


1.3 (2itl)(2J+j^) . . 

2.4 ‘(2i + 2)(2t4-4) 



it being remembered that for i = 0, the factor outside the s(juare brackets 
is 2. In the case of Venus and Mars is not far from 5 and the convergence 
is very slow. The plan adopted consisted in finding ea(;h coefficient to seven 
significant figures with seven place tables, so that the results might b() trusted 
to six significant figures : this involved the calculation of the portion inside 
the bracket to seven places of decimals f. 

A more rapidly convergent series can be obtained, especially for large 
values of i Expand in powers of a’, 

1 I 1 

and put = = 

1.3...(2i-l) 1 

2.4... 2i ’’21 + 2 


Then an easy transformation gives 

(23) A* = 2a' [(2i + 2) qi+iao - qi^iai - qi+^a^ -...]. 

A table of the values of a' was formed and also one of the logarithms 
of the qi\ the latter are rapidly obtaine<l by making a table for (2i~ l)/2i 
and then one for (2i + 2)qi and finally one for qi. In this way any coefficient 
A* can be easily and .securely found. 


Tisserand, loc. cit. 

t The phrase ‘ degree of accuracy ’ used in this Bcction refers to the number of significant 
figures in the results and not to the number of places of decimals. 
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We have the relation* 




2i- 

2z- 


2 

1 


eA^-^ - 


2t-3 

2{-l 


At-2^ 



which is ordinarily used to find A^ when A^~^ are known. But it is 
better to start from (23) with the highest value of i required and use the 
relation in the form 


(24) 




2i + 2 
27+1 




2i + 3 

27Ti 




The loss of accuracy is much less, and it diminishes with i\ in the case of 
V(jnus the loss amounts to about one significant figure after six successive 
applications. For Venus, the values of A^ were found from (23) for i = 43, 
42. 41 ; 30, 29, 28 ; 20, 19, 18; 10, 9, 8 ; 2, 1, 0; each triplet was tested by 
(24) and the intermediate values found from the same formula. An in- 
dependent test is obtained by forming the successive first, second, . . . differences 
of logA‘ which rapidly become constant when i is greater than about ten. 
This fact was used to correct tlie intermediate values where necessary. For 
Mars, a is a little smaller and the loss of accuracy less; but it was found 
more rapidf to calculate for t = 30, 28, 20 ,... 2, 0 from (23) and to find the 
odd values from (24) with i an even number (since e > 1), testing with i odd. 
For the other planets it was sufficient to find the highest values from (22) 
and follow down with (24) to i = 0, which was computed from (22) as a test. 

For the we have (Tisserand, loc, cit) 

( i ^ 7^>i) == « (^i _ (27 + 1). 

(2o) I ^ ^ 

1 i (B< - «'+•) = 2 (1 + + 1). 

from which any pair of values of the IP may be computed ; there is however 
a loss of accuracy. But we have (Tisserand, loc. cit) a formula which may 
be written 


(20) y7‘-*-2/A* + i7'+‘ 

and which enables us to calculate the 7?* with great rapidity, without any 
loss as soon as the two highest values of have been found ; moreover the 
loss in the latter disappears after the first step and the succeeding B'~ 
have the same degree of accuracy as the corresponding At For security 
the were calculated from (25) after each ten steps, and several were 
tested by the formula (Tisserand, loc, cit,) 




27 ~ 2 
2t - 3 


eB^-' - 


27 - 1 
27-3 


Bi-\ 


* See Tisserand, .Wc. CVl, , Vol. i. Chap. 17. Tisserand gives the formnlae for 
A\ 

t The calculations for Venus had been completed before Mars was undertaken. 
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The greatest advantage of (20) however is the rapidity with which the 
calculations are performed by means of it, for it does not re<juire the use 
of logarithms. 


Special values of the remaining functions O’, D\.,. were obtained from 
the formula (to be proved later on) 


(27) 


*+* ■^6'>S(S+ 1)^*’ 




which involves no loss of accuracy, combined with 

(28) ^ ■ 
for the intermediate values. 




+ a, 


(It is interesting to notice that (27) gives 

P! = P-i’ 


but no use was found for this formula.) 

The have now to be obtained. There are two problems to be 

considered : the first is where a table is to be made for all values of i up 
to some definite place (30 in the case of Venus* and Mars, 8 for Mercury 
and 6 for Jupiter). Wc require here a formula for special values of i, and 
when the coefficients for the two highest values have been found, the 
remainder are obtained from the obvious generalization of (28): 


(29) 


yTji— 1 


A./ 


+ /9 


i'i'l 


This method does not give the values for .? = 4 ; the formula for special 
values of i is required throughout. The .second problem arises when so few 
values of i are needed that it is not worth while to make a full table. Hence 
we require to develope formulae from which the can be obtained for 
special values of i, it being understood that the tables are always formed with 
the assistance of (29). 

In the following, the index is generally omitted, since it is the same 
in all the formulae now to be obbiined. 


Remembering that 

Ao“* = {1 ~ 2a cos (2^ — P) 4- a*} "2 =s cos i (T — P), 

we easily obtain, if for a moment we put 

Db, = s ( I - a*) bg^i - sbg, 

D^bg = + i^bg, 

{D-\-28fbg^4^bg+,-\-i'bg\ 


Except for A, B, which were required to f = 43. 
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J {D — 5 4- i)*^« = "h 

(D + * + i)«A=vA+. + i‘A, 

OL 

IJ^l3. = + 2 (s - i) D/S. + \i‘‘ - (s - /3., 

(D + 1 ffi, = 2s {se + ie') /8,+, + (i’ - + i) A- 


[sect. 


By definition, 


^«,r = 




jp ■ \daj 


A 


1 


= iD(D-l)(D-2)...(D-^ + l)^., 

whence, if k be any number independent of a, 

(.S3) (D + k) /3,,j, = (p + 1) A, P+, + (p + A:) 

(34) (D + A;)’^,,p = (p + l)(p + 2)/3,,p+, 

+ (p + l)(2p + 2A:+ l)ft,p+, + (p + fc)'‘^,,p. 

Again, from (30), 

(35) (D-S+ 4)“^.,p = D(D - 1) ... (D - p + l)(4s»a/8.+, + 1“/3.) 

— 46 '* ot 4 - 

Combining this with (34), after putting A; = — « 4- we find 
(30) p {p 4- 1) = 46“a 4- — p {2p — 2s) 

4 - [p -(p-S- J)=^) ffg,p-r. 

The material contained in these equations is sufficient to find all the 
^t,p rapidly, when the fig are known. 

We have, in fact, 

Ai = Je J5o — Ji4o» 

(;, = §€'i)o~iCo, 


2^3= aBo— il, 4-i’^o. 

22?3 ?= OolCq 4" Bi 4" (i^ — l)-^o> 
2O3 = 25aD, 4- 30, + (i* - 4) Oo, 


6 ^, = a (B, 4 2 ?o) - 6^3 + - 1 ) il„ 

GZl, = 9a (0, 4- Oo) - 2B, 4- i^B, , 

GO, = 25a (Z), 4- A) + 2O3 4- (i* - 1) 0„ 


12.4* = a ( A + A) ~ 15^3 4* (i* - 4) At, 
12A = 9a(03 4-0,)- 9A4-(»®-1)A/ 
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The required functions >0) are obtained from these without loss 

of accuracy. 

The functions {(D + 1)® — i''*} (D^ — are required to find the 

Pi, Pe, p 7 ,...(page 24). If these were obtained from (34) with the values 
of just found they would in general entail a loss of accuracy of the order 
of i units in the last calculated figure. 

Equations (31), (32) enable us to calculate them without loss when 
p = 0. When > 0, we have from (34), (36), 

[{D 4* ly — i‘\ == 45^a(^,+,,p 4- ^«4-,,p_i) 4 (/) 4- l)(2s 4- l)^,,p-|.i 

+ (5 + i)(2p~«4-f)ft.j)* 

In particular, 

j(i) + 1)0 _ Uo = + 3a) 0, - 2ii,„ 

{(D + ly - Co = -I- 10a) D. - 6C„ 

{(/) + 1 y - i‘] A,,= a(B^ + 7},,-,) + 2 (j3 + 1 ) + (2ji + 1 ) vlp, 

{{]) + 1 )• — i-1 Bp = 9a (Cp + C,,_,) + 4 (p + 1) Wp+i + 2 . ipBp, 

[(/) + ly _ is) Cp = 25a {Dp + + 6 (p + 1) Cp+. + 3 (2^) - 1) Gp, 

which are rapidly obtained, since the multiplication by a will have been 

made in finding the 

In a similar way we obtain from (31), (34), (36) for outer planets, 
(D“-i»)Ao= olBo, 

— i^) Bq = 9a(7o 4- 2Pi — • Bo, 

{D^ - i^) Go = 25aPo 4- 4(7i - 4>Co, 


(D^ — i^) Ap = a {Bp 4- Pp-i)> 

{D^ — i*) Bp = 9a {Cp 4- Gp—i) 4- 2 (j? 4 1) 4* (2p — 1) Bp, 

(P= - i^) Cp = 25a {Dp 4 Dp.,) 4 4 (p 4 1) 4 2 {2p - 2) 

It has been found convenient to alter certain of the formulae (21) a little. 
We have iA^ = iP""’ - 

A.pplying the first of these to the definition of p*, the second to that for 
L\ and both to that for G\ it is easy to show that 

( 37 ) 

^ 3 

BB. ^ 
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These, equally with (21), can be operated on by any function of D and 
consequently they are true when the suffix p has been attached to the 
symbols. 


Solution of a Dijfficulty. 


From the series for Apy or otherwise, it is obvious that — is of 
the same order of magnitude in general as A^. Hence when ^ is large and 
this function is calculated from the tables it appears as the difference of two 
nearly equal numbers, and a considerable loss of accuracy results. The same 
fact is true of all the functions Apy Bp , ..., but for a given value of i the loss 
of accuracy is smaller as p increases and as we go along the series A,B, 0 , — 
Although the number of places used in the tables has been found to be 
sufficient for the purposes of this essay, I shall give a method showing how 
the difficulty may be overcome. 


Wo have (Tisserand, loc, cit.) : 

. . 2 . sin-' ylr 

A^ = aM - .1 

TT •'o (1 — a^sin-‘v|r)2 

Hence ' (D — = - I " r - 

TT Jo (I — a'sin^\!r)--i 


(I) -i- l)A'+> = -- a‘+» r dyir. 

TT Jo (1 — sin® 


Therefore 


I (1 - a® sin® 
(D - i) A^ = 4- a (Z) - t - 1) 


= a A 4- a®A^+® 4- ... 4- 4- aJ (D - i - j) 

The values of A\ 7)A' = A^^ have been calculated to seven significant 
figures, and the ratio of (D — i — to (D — i)A^ approaches the limit 
a as i increases. Hence if p is the number of units loss of accuracy in the 
last significant figure of {D — i —j) A^^J, the number of units loss in (D — i)A^ 
is pa^ approximately, which can always be reduced to unity or to a gain by 
taking j large enough. For values of i less than i' we must take a-^ < t if we 
are to lose no accuracy in {D — i) A\ In the case of Venus, with i' = 30, this 
gives j = 5, i = 25. 

It is obvious that this principle can be extended to any of the functions 
previously used by applying the preceding formulae, but since the practical 
applications did not demand it I shall not develope the formulae here. 

It does not seem possible to obtain a general formula for {D — i) 
in terms of p* without loss of accuracy otherwise than by the principle 
of the method outlined here ; i.e., by referring forward to some higher value 
of i. 

The corresponding difficulty in Radau’s work appears in the form 
A' - olA'-\ a* - 2aA*-' 4- o®A'-*, .... 
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Calculation of the Goejfficients Mi and of their Derivatives with 
respect to n. 


The most arduous part of the work has consisted in the computation of 
the derivatives of oc^ — + ^ixy = Uo“ with respect to a, to the degree of 

accuracy which was demanded. The only complete literal exj)an8ion in 
powers of m is that given by Delaunay, and his results suffer from two serious 
defects for this purpose, the one in the fact that the coefficients for the 
parallax on which the desired quantities chieHy depend are not taken far 
enough, and the other in the doubts as to the degree of convergence of 
the series along powers of m. Both may to a certain extent be remedied 
by a judicious use of the numerical results given in my Lunar Theory. So 
far as I know, the only method for calculating these derivatives from 
numerical results in which the value of m has been substituted, without 
having recourse to algebraical series in powers of niy is one which I gave in 
1903*; this method, although tedious and difficult for numerical application, 
has perforce been adopted in the absence of any other plan. 


The method essentially consists in finding three independent relations 
between the coordinates and their derivatives with respect to the six elements, 
expressed in the form of algebraic or differential o(|uations. From these 
three relations it is possible to obtain the derivatives of the coordinates with 
respect to one of the elements in terms of the others. No approximation 
processes are necessary — only multiplications of series and a cjuadrature for 
each coordinate. Further, the results naturally appear in the form in which 
they are required in the disturbing function, namely, as derivatives of the 


rfc, 


rectangular coordinates, and the value of 


, on which the principal term in 

the equations of variations mainly depends, also arises naturally. The order 
with respect to e, k to which the results can bo calculated is one less than that 
to which the coordinates have been obtained, but it was not found necessary 
to proceed to high orders. 


The formulae are constructed to solve the problem : Given the lunar 
coordinates and their derivatives with respect to Ca, c.,, Wiy Wa, to find the 
derivatives with respect to n, it being supposed that the coordinates are 
expressed in terms of these six elements. My lunar theory enables one to 
get all the data of the problem accurately. It is expressed in terms of the 
constants e, k, which are left arbitrary; but the second and third of 
equations (8) together with the values of c,, Cj in terms of these constants 
enable one easily to find the derivatives with respect to c,, c» from those with 
respect to e, k, or vice versd. The method given for performing the calcu- 
lations was adopted almost without change. The results will be found in the 


Tram. Amer. Math. Soe. Vol. iv. pp, 234 — 248. 
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collection of tables. The details by which the functions chiefly needed, 
namely, j j 

were found directly instead of the functions for which the formulae were given, 
namely, 

a / V O' / \ 

dnS ^' 

will be omitted*. 

The greater part of the work consisted in multiplications of series of the 
same nature as those necessary in my lunar theory. The calculation of 

Uo", \x,z 

comsists entirely of such multiplications. The multiplications by 
a'- 

rp'i * ^ 


are (piickly performed. Certain coefficients of orders higher than those given 
in the tables were re(|uired. These were obtained from the complete results 
already ((uoted for u,,, z'. their derivatives with respect to n were multiplied 
by such small coefficients that a calculation of them was obviously un- 
necessary. Thus in all this work the efiect of the solar perturbations in the 
moon 8 motion on the coefficients of the planetary terms has been fully taken 
into account. It may be stated here that the derivatives with respect to n 
are only recpiired either in the terms of comparatively short period (say, a few 
years or less), or in those terms which do not contain the argument Wi, 
A glance at the final form of the eipiations of variations will make this state- 
ment evident. 


* Some orrora in the formulae of the paiHjr referred to were found. The integrated equation 
(C) p. 246, requires an added constant, contrary to the statement there made : this constant may 
1)6 determined by the equation numbered (14) in the pai)er ; and the sign before the last term of 
(/•;) p. 247 should be changed, as well as the first negative sign in equation (15), p. 244. 
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A SIEVE FOR THE REJECTION OF INSENSIBLE COEFFICIENTS. 


The periodic inequalities due to the direct action of the planets arc 
usually divided into two classes, which have received the names ‘ long-period ' 
and ‘ short-period.’ These names are somewhat misleading in this connection, 
because there is no sharp line of division. I shall divide them into primary 
and secondary inequalities : the former are those which arise from the substi- 
tution of Wx for the non-periodic term of V ; the secondary inequalities 
will be defined as those which arise from the substitution of the variable 
values of the elements in the periodic terms. 

The majority of the primary inetjualities are of long period — a year or 
more ; but there are classes of them having sensible coefficients with periods 
of a month or less. Nearly all the secondary inequalities are of short 
period. Thus the primary inequalities in the longitude are obtained fronj 


and the secondary inequalities from 




1 


0^2+ , - ow., 

/ dWi^ 


^dV . ^dV . 

-f- , - or? -f J — 0C3 -f- j — OC‘3, 
dn dc.^ f/Ca 


The coefficients of the secondary inecjualities will be smaller than the 
coefficient of the corresponding primary ine(|uality, in general ; an examination 
of the final form of the equations of variations and of the values of the co- 
efficients in the derivatives of V shows this immediately. An exception will 
occasionally occur in the terms 


dy 

dc<i 


Bc-if 


dV 

dw,. 


hw.. 


The largest term in 
hct is 


dV 

dc^ 


(see page 12) has a coefficient 40 ; the coefficient in 
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that in Bwi is 


C 




dV 


Hence the ratio of the largest coefficient in Sci to that of Sw, is 


since 


1 2S si., 1 

~ 2 10” • X, • , ,/'s A, 
^/X. • A 


2 

/ 10 ' ■ 


14l2 


10“^J+ 2^' 
H A 


A ^1* 

As IS 
A 


never very small (it is in most cases a little greater than unity) 


the value of p will be large only if the two terms of the denominator are 
approximately e(|ual in magnitude and opposite in sign. The sieve which 
follows does not contemplate this ‘ accidental ’ approximate vanishing of the 
denominator, so that it would not reject the inequality in for this cause. 
Hence the maximum value of p would appear from the sieve as the greater 
of the two values 



144.9 s 
10”X, 100 ■ 


since X, =4 4- ()148G4— *0087444 and its least value when 4^0 is taken to 
be *01486 in the sieve. The terms in which the former is the greater are all 
calculattjd without reference to the sieve. Those in which the latter is the 
greater demand that s > 100" for the secondary inecpiality to be rejected by 
the si(^ve for when it ought to be retained. The sieve shows that such 
inequalities can only occur with the lunar arguments 2D — 24 D — 4 and 

dV 

these have been separately examined. The term is of the same order 

approximately, and the same argument applies. 


A similar argument applied to the terms containing 21) ~ 2F, D — F shows 
that the sieve would not rtqect the secondary inequalities in longitude as 
long as s<2000", and as the upper limit for the sieve is s=3500", a secondary 
ineipnility with a maximum coefficient of 0"*02 might have been rejected by 
it. But the ineipialities for 5 > 2000" were separately examined. 

The secondary inequalities in the latitude of terms containing 2D — 2F, 
D — F are more important than those in the longitude. A similar 
argument applied to them shows that the primary might be neglected by 
the sieve if s > 60" when the secondary should be retained. Consequently 
a separate examination of these terms was also undertaken. 

We have then to examine the primary inequalities. Those which have 
periods of less than a year do not need a sieve. The terms in the moon’s 
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coordinates which, in combination with the planet’s coordinates can give 
sensible terms, are few in number ainl they are taken one by one according 
to the sizes of their coefficients with all possible terms arising from the 
planet s coordinates until they become insensible : the number of the moon 
terms thus required is quite small, as the results show, and any omissions 
would arise only from numerical errors. The sieve whicli follows is therefore 
constructed for periods greater than a year. Some limitation had to bo 
placed on the maximum period to be included ; it was taken at 3500 years, 
corresponding to a value of .v equal to 1". But values of s less than I'" which 
arose from the equations constituting the sieve were also included. 

The results show that though th(‘xe will be inequalities of periods greater 
than 3500 years, their order is so high that it is extremely doubtful if their 
coeOicients would be sensible. In any case the period is so long that they 
would scarcely be observable within historic times. It is to be noted that for 
very long periods it is better to expand such terms in powers of t when 
we wish to examine their effect on the observed position of the moon. 
Thus : 

Osin (qt + q') = cto + otit-{’0iit^+ .... 

The first two terms would then constitute an addition to the expression for 
the moon’s mean longitude, which, being one of the observed quantities, would 
not be affected : the values of «<,» «i could hardly be large enough for this 
change in the meaning of to affect the periodic terms. The succeeding 
terms would constitute an addition to the secular acceleration, but it is quite 
improbable that there is any coefficient great enough to affect it by so much 
as 0"T per century. 

It has not hitherto been found possible to construct any kind of analysis 
which will give a criterion for the magnitude of th(^ coefficients of the 
planetary terms without having recourse early to numerical developments. 
There are many thousands of ine(|ualities who.se periods suggest that their 
coefficients might be sensible but which on computation are fotnid to be 
exceedingly small. Even the roughest approximation would be a laborious 
process, and it seemed advisable at the outset to construct formulae which 
could be rapidly applied. Fortunately, for the two most troublesome planets 
— Venus and Mars — numerical developments were available for this purpose. 

Newcomb has given • the expansions of in cosines and sines of 

iT—jF up to 1 = 29, with a sufficient number of values of 
numerical values of all quantities have been substituted at the outset, so that 

• Amer. Eph, Paper$t Vol v. pt. 3. 

t The order of any term with reference to the eccentricitiea of the earth and planet and the 
inclination of the planet is 1 1 j (. 
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there is no possibility of using the method of Section ii. But as the argu- 
ments are developed in multiples of the mean longitudes it is still possible to 
differentiate or integrate with respect to these, on the understanding that we 
consider the arguments as functions of the mean longitudes and of the longi- 
tudes of the perigees and of the node (not of the mean longitudes and mean 
anomalies as in Section II.). For Mars the development had not been carried 

beyond i = 1 1 for , and not completely to that point ; but a process of 

extrapolation and a knowledge of the general law of decrease enabled one to 
extend the table for the purposes of a sieve sufficiently far : a liberal margin 
was left for errors in this extrapolation. For the other planets, the general 
law of decrease, which approximately holds when i is sufficiently great, namely 
= a, is so soon reached owing to the smallness of a, that there was no 
difficulty in forming a rough table. 

There are two problems to be dealt with : first to form a table from which 
all possible periods can be rapidly obtained ; second, to obtain formulae for 
finding easily the order ot magnitude of the coefficient which will result 
in the moons longitude. The process which follows is divided into several 
parts : — 

(а) A formula for the coefficient in longitude when we know A 
(the coefficient from Newcomb’s tables), s (the number of seconds in the 
daily motion of its argument), k (the order of the term with respect to the 
lunar eccentricity and inclination). Modification for the planets other than 
Venus. 

(б) A formula for the order with respect to e, y associated with a given 

multiple ot tei, in order to obtain k, 

(c) A method of constructing the period corre.sponding to each multiple 
of k to be considered, to obtain s. 

(d) The terms which were not rejected. These will be found with the 
numericiil results. 

The only part ot the actual computations which required special care was 
that giving the periods. It was found, however, throughout the numerical 
work that, the plan once formed and in working order, the size of the 
coefficients could be almost guessed at a glance, so that the computations 
were themselves of the nature of a test of accuracy, and errors tended rather 
to make the coefficient too large than too small. The test was gone through 
twice. The detidls which resulted in the table of coefficients retained by the 
sieve for accurate computation will be omitted. 
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Comtructimi of the Sieve. 

We have, as in Section ii., 

( 1 - 7 "^) r" cos ( V" ~ F) + f'^r" cos ( V' + V' - 2h"), 
7 ; = (1 — 7 ^ 3 ) F sin ( F" - F') — y'^r" sin ( F' + V' — 2 /i"), 
f = 27 " sin ( F' ~ /i '), 

A* = F-^ - ~ 2r'f. 


Hence, putting 


we obtain 


,, rfr ,,, , dV' 

A^J^, = -^E-r,N, 


n_N+JE 

^«+2 


1 /^A.\ ; 

s dT Uv 


Multiplying the latter equation by l^E—r)N, f 4 - successively, and 
making use of the equation A'^ == {^ 4 - F + form* 


we find, after a few simple operations, 


f- - f ^ N'^ - E'^ /I _ __ 2 (7;jV - f A’) d / 1 \ 

' “ jV‘i 4. A'i \A-'' A'‘^V .9 (jV* + 7 ^^“) V^V 

2 iVA’ / 1 _ mN-^vE) d i\\ 

A''^ jv-i q: ( jV2 ») ■ 

Now the principal terms in the expansion of N, E in powers of e' are a\ 
ae\ and the equation (38) shows that when the multiple of T considered in 

the expansion of -- in cosines is not very small, is small compared with 
positive). Hence we shall evidently be taking the worst case if we 


consider 


r-F 2|5 


of the onlers 


^ d (\\ 

3 c^rUv’ 


4f d /1\ 
■3 dT Uv ‘ 


Further, 


d / 1 \ (7 / 1 \ 1 rff 
AV" 7 'rUv l^^ dT 


= order of 4, 


’ aJ’ 


* The symbol i is only tu»ed In this connection on this page, and it will not be confnaed with 
the « in the eqnations of variations. 
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which, for the worst case — that of Venus where a'- = 2 a^ approximately — 
gives 

The order of tj is approximately the same. 

Consider next the function 

I 

A* A® ' 

Here f is of the order 2 y", Suppose that these functions be expanded in 
powers of and then in cosines. For Venus 7"' < and therefore for a 
given multiple p of T in which p < P\y the terms factored by 7"* will be less 
than those independent of the value of j>, which would make them 
greater, is seen in the course of working the method to be too large to need 
consideration. Hence for coefficients which do not contain 7"^ as a factor*, 

it is sufficient to take : for those which contain y'- as a factor we should 

take , but this being less than within the limits, the latter was 
used. 


The first term of the disturbing function R in Section ii., when written in 
terms of the real variables, is 


R = m" |i (?•“ - 3^0 ) J - y') (' 


r-v 

A» 




in which the first factor is the part depending on the moon. If we take 
^ = 7 = tV instead of their actual values general take no 

account of the fact that wt* shall not be making the coefficients too 

small if we consider — 2/'^ 2.ry, 2 zXy 2 zij as of the order 

where k is the order of the term considered with reference to e, 7. Remem- 
bering that a combination of two sines or cosines introduces the factor we 
find the following results, in which the factors have been increased slightly to 
make them suitable whole numbers. 


First term of Ry order - ^ 10"*^ when 7"* is not a factor, 

0 A’ ’ ’ 


A» 


10 ^ when 7"“ is a factor. 


Second and third terms of R, order ^ jL 10“^ 

5 al A* 

fourth term of R, order ^ ^ I0~*. 

10 A® 


The parte depending on 7"^ always diminish the coefficient. 
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Next, the primary inequalities are found with 

+{; J (4a»)] sin {qt + 5'). 

where 

1 m" 1 in" 

~ 4 ~ 4 ^ cos + q't + 9'"), 

\ = -‘ h 4- 0148612 - 0037441,, 

There are several cases to be considered. 


Case 1. 1*1 = i, = i;, = 0. The second term of S?e, with the first term 
of R since the other terms of R will have the factor at least; 

I w. I . , , , , 


^{An^)\^^\A\, A; = 0. 


Coefficient in is of order 


• (f'^\ r 


coef. in 


^ not a factor, 


7"' a factor. 


Case 2. ii — 0, f, even. The first term of combined with the second 
or third terms of R ; 

Coefficient in 8?e, is of order 


tn" 


(coef.iu-4«;). 


Case 3. u = -h, i, ^ 3, i, even. The first term of combined with 
the first term of R ; 

Coefficient in Swi is of order 

1012.8-i: 

The other cases, depending on the way in which Wj, w^, w, enter into the 
argument, are treated in like manner. 

I give next the final form of the equations used in the case of Venus. 

Here Newcomb tabulates 24^,, 24^^. If are 

m 4(Jo,Ui/U ^ ^ 

the coefficients derived from Newcomb’.s tables for these functions, p the 

multiple of T present, C the coefficient in Bwi expressed in seconds of arc, 

we find the following expressions for logjoO; the arguments are put into 

Delaunay’s notation {il + ig -f i'Ci). Additions have been made to the 

numbers for simplicity in the results ami for parts arising from the terms 

in the moon coordinates containing e\ 
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Argument 
/, h absent 
I absent, i + i'\ even 

„ i + i", odd 

I present, i = - 0, 1 

„ i' even 
„ mid 


LogC 

1 *5 — log 8 + log jVg, 

greater j4*5 ~k - 2 log 8 + log p + log N’s, 
of l3*5 - A: - 2 log s + log iVs , 

3*5 — ^ — 2 log « + log iVs, 

Q'O — k — 2 log 8 + log 

5-0 - A; — 2 log 5 + log/? + log iTg, 

5*0 — A; — 2 log 8 + log iVj . 


From the second term of R, which produces terms containing the factor 


a 

a'’ 


by a similar procedure I find 


Argument 
1 = 1 , 

1 = 3. 


% + X even 


' + i" mid 


LogC 

2*3 — k—2 log 8 + log iVg, 

2*0 — A; — 2 log « + 2 log p + log 
1 *5 — A: - 2 log 8 + log p + log N^. 


tor the other planets, similar expressions can be obtained. Except in the 
case of Mars as far as 1 = 11, no tables giving iV'g, are available. It is 

therefore necessary to construct a method for the coefficients in “ . . As 

A'’’ 

stated above, the parts of these functions independent of the eccentricities of 
the earth and planet can be expanded rapidly for a rough approximation. 
Ihe terms which contain the eccentricities jis factors are in general of the 
orders {pe or {pe^i (where p is the multiple of T present, q the order with 
respect to e, e") compared with the order of the term with argument 
p {T - P), as a ref(*rence to Leverrier’s development shows. A margin was, 

however, left for the occurrence of numerical factors independent of — , p, 

but depending to a small extent on q, Jupiter was the only planet in which 
the approximations to the coefficients came close to their correct values ; and 
this occurred where the approximation was comparatively simple, namely, for 
small values of p : a smaller margin has been necessary. 

The coefficients which the sieve gave as greater than 0" 01 will be found 
in the collection of numerical results. 


With a knowledge of the greatest values of iVg, N, (which arise with the 
lowest values of p) it is a simple matter to find the maximum value of s to 
be considered in any given ci\se. In fact, the maxima of s are found by 
eejuating the above expressions to -2, since the least value of C to be 
considered has been taken at 0"*01. Many limitations which it is not 
necessary to specify in detail appeared during the progress of the work ; e.g., 

if the multiple of I present in the part of R independent of — , is 2, it is only 

necessaiy to consider values of p lying between 25 and 50, and then only 
those in which the sum of the multiples of P, - T is less than 5* 

• These statements are always made with coordinates referred to axes following the mean or 
true place of Uie sun. 
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Lowest orders with reference to e, h associated with the Lunar 
A7'gu7tie*fts\ 

If any argument contains w,, w^ in the form* 
iitOi 4- 4- 4^8 = il + i'g 4- t'h, 

so that 

I “ ^'1 > t ~ 4 " 7 ^ 2 1 1 “ ^>1 4 " ^8 4 " ? 8 * 

then it is easily seen from the known properties of the solution of the problem 
of the moon’s motion, that the lowest orders are as follows : — 

(40) 1 1 I 4- 1 i' — i" I for i, i" even, 

(41) 2 4- h — I 4- 1 1' — i" I for i\ i" odd, 

(42) 2 4- 1 1 — i' I 4- smaller of 1 1 ' — i" 4 1 1 for i' even, i" odd and positivef, 

(43) 4 4- 1 1 1 4- smaller of 1 1 — i" 4 1 | for i' odd, i' even and positivef . 

The first expression (40) represents the orders for all terms in longitude 

and parallax or in oc, y which are independent of ; (41) for terms con- 
taining odd powers of as a factor; (42) for terms in latitude or in z 

independent of , and (43) for those containing an oild power of as a 

factor. Unity has been added in the expressions (42), (43) since an odd 
power of 7 in R is always associated with an odd power of 7 ", and 
conversely. 


Method for finding the Pemods to be Examined, 

The lunar arguments r/, h are first considered. These are divided into 
classes according to the multiple of I present ; when the multiple is greater 
than 3 , the value of p is so large that the coefficients are (piickly seen 
to be insensible ; in fact, there are none retained by the sieve with the 
multiple 3. 

Consider any multiple i'g 4 i"h of g, h. All the facts needed concerning 
a given multiple of g, h can be inserted on a sheet of squared paper by taking 
i’, i" as the coordinates and inserting in the scpiares: ( 1 ) the number of 
seconds in the argument ; ( 2 ) the order with respect to e, 7 of the coefficient 
as derived from the expressions (40), (41), (42), (43) associated with the 
argument ; (3) the multiple of T associated with it in the greatest term 
of this order in the moon portions of R\ that is, the term whose coefficient 
in these functions does not contain e' as a factor. The specimen given 

* The results were worked out with I, g, h instead of with ic, , ii>„ ir,. 
t The positive sign is attached for oonvenienee: the sign of the coefficient can always be 
arranged so that V’ is positive. 
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on page 47 includes all terms of the seventh order. It is evident that only 
two quadrants are required. 

Similar charts were constructed for + g ’¥i"h, ... but 

that it was not necessary to insert the values of a will presently appear. 

For the planet arguments the table which follows was formed to give the 
periods of pT — up to the first very long one, excluding periods greater 
than that of T. The last three* given are the arguments of very long 
periods : it is not necessary to go further, since all the periods greater than a 
year are either multiples of those given or are combinations of the last three 
with the others ; higher multiples of T, V make the difference 1 p — | 

(which gives the order of the coefficient with respect to e\ e", 7 ") so great 
that the terms would obviously become insensible if 5 > 1 We thus have 
all the periods of arguments independent of I, g, h. This particular table 
also illustrates the plan on which the sieve was worked. 


Arg. 

s 

\N\ 

logic] 

For Cal. 

T 

+ 3548 

110 

- -9 -'5 

X 

12r-8F 

- 3563 

•03 

-3-5 


47^-3r 

-3110 

14- 

- -9 --5 

X 


+ 3096 

•04 

no 


or-er 

-2672 

•3 

-24 


iT-2V 

+ 2657 

2-4 

-1-5 

X 

r- V 

- 2220 

230- 

+ -6 -*5 

X 

12r-7F 

+ 2204 

•002 

no 


i)T-iV 

-1782 

2*3 

-1-3 

X 

7T-iV 

+ 1766 

•3 

-2-2 


IIT-TV 

-1344 

•03 

no 


2r- V 

+ 1328 

14- 

- -5 --5 

X 

3T-2V 

- 891 

13- 

- -3 --5 

X 

10T-6V 

+ 876 

•03 

- 2-9 


8T^5V 

- 454 

•25 

- 1-5 

X 

5T--3F 

+ 438 

2-3 

- -7 

X 

13T-8F 

- 15 

•004 

-2-0 

X 

26T-16F 

- 30 


no 


39T-2iF 

- 45 


no 



Formula : log 0 = 1’5 — log s + log N^. 



Max. of 9 

200 

80000 

20 

8000 

2 

800 

•2 

80 

•02 

8 

•002 

1 


— *5 to be added to log | C j when the multiples of 2\ - V differ by less than 2. 

* These three are multiples of the first of them, but it is oonvenient to have them written 
down. 




Multiples of h 


TABLE FOR THE PERIODS AND ORDERS 


Multiples of g 

1 2 3 4 5 6 
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For arguments containing g, h but not ly we proceed according to the 
order given on the preceding page (for h\ beginning with the lowest, which 
is 2. There are three terms of this order with arguments : 

2^ -j" 2^-, 2/t, hy 

and motions 

801", - 382", - 191". 

Taking the first, which arises principally from the terms 2 x 1 / in R 

in the form 2D — 21, we examine by means of the second formula on 
page 44 all the terms which have values of s less than 3500, obtained by 
adding ± 801" to the values of s in the table given on page 46. For 
example, the argument 

2g + 2h — lOr + 6 F gives s = ~ 74. 

To obtain we write it 

2D-2^-8T+6F, 

where D is Delaunay’s notation for the difference of the mean longitudes of 
the sun and moon. The value of obtained from the term in Newcomb’s 
table corresponding to — ST-i-GF, is extracted, and log (7 is then obtained. 
Before finding the terms of a given order, the maximum value of s needed for 
that order was found and the process was followed until the maximum value 
of 8 became less than 1". 

For terms containing the first multiple of I, the most convenient one of 
very long period is taken, say ^ + 16T ~ 18F in the case of Venus, for which 
s= — 13". This, after the maximum of s has been found, is combined with 
the motions in the table of arguments pT — giving a list of terms to be 
examined. The next first order term is 2D — I, which is treated in the same 
manner; the starting term is 2D — ^4*8T— 12F, for which s = — 87" The 
same process is followed with every moon argument : it proceeds according to 
the orders of the coefficients with respect to e, 7, and is continued until an 
order is reached for which the maximum of s is less than 1". 

The charts and tables which were constructed, together with the expres- 
sions constituting the sieve, permit one to review rapidly all the possible 
terms without any fear that some may be omitted. All the major planets 
from Mercury to Neptune were examined. For the minor planets a different 
method of treatment would be necessary: the principal effect is probably 
that of a circular ring of matter at the average mean distance : but the mass 
would be an exceedingly doubtful quantity. 

It should be added that the numbers given in this section are only rough 
and ready approximations. Accurate values were afterwards used in the 
actual computations of those terms retained by the sieve. 



SECTION V. 

AUXILIARY NUMERICAL TABLES. 


Epoch 18500. 


Arg. 

Daily motionn 
of arguments 


Longitudes at Epoch 

Perigee 

Node 

U)^ 

47434-891 

Mercury 

75^^ 07' 19" 

46° 33' 12" 

w.. 

400-923 

Venus 

129° 27' 34" 

75° 19' 47" 

w-i 

- 190-772 

Earth 

100“ 21' 40" 


Q 

14732-420 

Mars 

333° 17' 55" 

48° 24' 01" 

V 

5767-670 

Jupiter 

11° 54' 27" 

98° 55' 58" 

T 

3548-193 

Saturn 

90° 06' 40" 

112° 20' 51" 

M 

1886-518 




J 

299-129 




S 

120-455 






Eccentricity 

Inclination 

Sine half inclin. 

1 0) 

m' 

m" 

Moon e = ... 

•10955 


k -044780 



. „ e= ... 

-054906 


y -044887 



Earth 

-016772 





Mercury ... 

-205604 

7° 00' 07" 

•061066 

1-5878216 

6000000 

V enus 

•0068446 

3° 23' 35"-3 

•0296063 

1-8593374 

408000 

Mai-s 

•093261 

1° 51' 02" 

•016149 

•1828960 

3093500 

Jupiter 

•048254 

r 18' 42" 

•011466 

•7162374 

1047-35 

Saturn 

•056061 

2° 29' 39" 

•022 

•9794957 

3501-6 


a difF. of masses of Earth and Moon 

a, = — . — *0025000. 

' a sum „ „ „ 


BFL 


4 
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a'Ap\ Ventis. 


i 

p = 0 


p = 2 

p = ‘S 

P = 4 

0 

2-386373 

1 -18898 

2-01958 

3-56717 

7-14018 

1 

•042412 

1-64375 

1-97018 

3-61813 

7 15763 

2 

•527578 

1-49718 

2-23897 

3-67099 

7-24792 

3 

•323341 

1-25777 

2-38378 

3-90485 

7-38420 

4 

•206787 

1-01817 

2-37596 

4-18715 

7-66163 

n 

•135585 

•806423 

2-25093 

4-40529 

8-06778 

6 

•0903733 

•629509 

2-05234 

4-50184 

8-.52231 

7 

•0609449 

•486319 

1-81709 

4-46256 

8-9;U31 

s 

•0414597 

•372759 

1-57213 

4-30024 

9*21822 

9 

•0283961 

•283965 

1-33525 

4-04108 

9-33621 

10 

•0195545 

•215254 

1-11704 

3-71519 

9-26837 

11 

•0135258 

•162505 

•922804 

3-35128 

9-02179 

12 

•00939037 

•122264 

•754294 

2-97337 

8-62021 

13 

•00653977 

•0917207 

•610989 

2-60025 

8-09624 

14 

•00456677 

•0686356 

•491048 

2-24529 

7-48568 

15 

•00319644 

•0512491 

•391967 

1-91729 

6-82294 

16 

•00224186 

•0381937 

•311003 

1-62103 

6-13881 

17 

•00157519 

•0284159 

•245-153 

1-35848 

5-45869 

18 

•00110854 

•0211093 

•192801 

1-12945 

4-80230 

19 

•000781251 

•0156602 

•150799 

•932354 

4-18393 

20 

•000551305 

. -0116035 

•117495 

•764686 

3-61295 

21 

•000389493 

•00858810 

•0912272 

•623494 

309465 

22 

•000275467 

*00634988 

•0706073 

•505651 

263101 

23 

•000195012 

f 

•00469065 

•0544896 

•408070 

2-22155 

24 

•000138178 

•00346202 

•0419391 

•3-27837 

1-86399 

25 

•0000979878 

•00255320 

•0322001 

•262283 

1-55486 

26 

•0000695400 

•00188159 

•0246666 

•209030 

1-28999 

27 

•0000493857 

•00138571 

•0188559 

•165996 

1-06486 

28 

•0000350954 

•00101988 

•0143859 

•131384 

•87493 

29 

•0000249553 

•000750191 

•0109555 

•103668 

•71570 

30 

•0000177550 

•000551514 

•00832894 

•0808952 

•58391 
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a'Bp\ Venus, 


i 

)i 1 

O ' 




p = 4 

0 

7-22787 

40-7577 

162-839 

574-438 

1889-97 

1 

6-41712 

39 9376 

162-009 

572-r)58 

1886-46 

2 

5-34305 

37-4703 

158-899 

567-202 

1875-65 

3 

4-30681 

33-9489 

153-053 

557-873 

1857-47 

4 

3-40299 

29-9236 

1 44 596 

r)43-770 

1831-35 

f) 

2-65251 

25-8036 

134-045 

524-376 

179fi-17 

6 

2-04714 

21-8593 

122-086 

499-716 

1750-67 

7 

1-56803 

18-2494 

109-417 

470-353 

1693-91 

8 

1-19391 

15-0509 

96-6472 

437-240 

1625-63 

0 

-904676 

12-2853 

84-2627 

401-549 

1546-42 

10 

-682782 

9-93951 

72-6128 

364-502 

1457-58 

11 

•513587 

7-98023 

61-9221 

327-246 

1361 05 

12 

-385215 

6-36440 

52-3111 

290-774 

1259-10 

13 

-288217 

5-04589 

43-8189 

255-885 

1154-16 

14 

-215180 

3-97965 

36-4253 

223-167 

1048-60 

15 

•160347 

3-12409 

30-0696 

193-016 

944-563 

16 

•119287 

2-41218 

24-6664 

165-649 

843-909 

17 

•0886081 

1-90189 

20-1175 

141-143 

748-121 

18 

•0657309 

1-47604 

16-3210 

119-461 

1 658-314 

19 

0487009 

1-14196 

13-1767 

100-483 

575-238 

20 

•0360433 

•880955 

10-5906 

84-0315 

499-324 

21 

•0266486 

•677817 

8-47687 

69-8952 

430-720 

22 

0196845 

•5-20254 

6-75904 

57-8446 

369-349 

23 

•0145280 

•398422 

5-37015 

47-6466 

314-956 

24 

•0107140 

•304484 

4-25252 

39-0734 

267-158 

25 

•00789548 

•232245 

3-35707 

31-9104 

225-485 

26 

•00591455 

•176824 

2-64251 

25-9593 

189-415 

27 

•00427940 

•134402 

2-07441 

21-0408 

158-405 

28 

•00314773 

•101996 

1-62429 

16-9955 

131-913 

29 

•00231406 

•0772887 

1-26880 

13-6833 

109-411 

30 

•00170032 

•0584848 

•988870 

10-9828 

90-4019 


4—2 
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[sect. 


a'Cy. Venus. 


i 

2) = 0 

i>=l 

p = 2 

0 

44-8769 

534-588 

3768-39 

1 

43-6375 

527-704 

3741-20 

2 

40-5989 

507-963 

3660-39 

3 

36-5135 

477-744 

3529-34 

4 

31-9852 

440-064 

3354-28 

5 

27-4388 

397-947 

3143-75 

6 

23-1435 

354-052 

2907-47 

7 

19-2502 

310-509 

2655-40 

8 

15-8260 

268-888 

2396-86 

9 

12-8837 

230-239 

2139-95 

10 

10-3980 

195-177 

1891-28 

11 

8-33082 

163 976 

1655-86 

12 

6-63168 

136-655 

1437-18 

13 

5-24910 

113-060 

1237-38 

14 

4-13379 

92-9238 

1057-42 

15 

3-24075 

75-9167 

897-406 

IG 

2-53032 

61-6831 

756-725 

17 

1-96836 

49-8668 ' 

634-298 

18 

1-52G09 

40-1283 

528-726 

19 

1-17959 

32-1543 

438-447 

20 

-909213 

25-6635 

361-821 

21 

-699011 

20-4082 

297-239 

22 

-536133 

16-1740 

•243-145 

23 

-410305 

12-7778 

198-107 

24 , 

•313370 

10-0649 

160-803 

25 

•238882 

7-90608 

130-066 

26 

•181778 

6-19418 

104-853 

27 

•138097 

4-84113 

84-2627 

28 

•104749 

3-77495 

67-5142 

29 

•0793390 

2-93721 

53-9426 

30 

•0600106 

2-28070 

42-9841 
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i 

€L 

31 

+ -0000126389 

32 

900151 

33 

641390 

34 

457215 

35 

326061 

36 

232619 

37 

166018 

38 

118527 

39 

846489 

40 

604735 

41 

432154 

42 

308914 

43 

220879 

44 

157975 

45 

113014 


i 


31 

+ -00172640 

32 

126734 

33 

929960 

34 

682120 

35 

500139 

36 

366576 

37 

268590 

38 

196732 

39 

144055 

40 

10.5451 

41 

771713 

42 

564604 

43 

412972 

44 

301987 

45 



Verms. 


i 


a'Dj 

0 

338-022 

6140-05 

1 

333-468 

6086-23 

’2 

320-567 

5928-97 

3 

300-989 

5679-86 

4 

276-750 

5355-68 

5 

249*813 

4975-76 

6 

221-872 

45.59-78 

7 

194-267 

4126-04 

8 

167-972 

3690-36 

9 

143-625 

3265-.59 

10 

121-595 

2861 -.50 

11 

102-033 

2484-89 

12 

84-9390 

2140-01 

13 

70-2012 

1828-94 

14 

57-6434 

1552-09 

15 

47*0520 

1308-57 

16 

38-1990 

1096-59 

17 

30-8580 

913-798 

18 

24-8142 

7.57-498 

19 

19-8702 

624-874 

20 

15-8493 

513-126 

21 

12-5965 

419-566 

22 

9-97758 

341-697 

23 

7-878,50 

277-235 

24 

6-20282 

224-141 

25 

4-87014 

180-612 

26 

3-81398 

145-080 

27 

2-97966 

116-193 

28 

2-32255 

92-796 

29 

1-80647 

73-913 

30 

1-40221 

58-724 




4785-68 
4743 14 
4618-95 
4422-57 
4167-46 
386907 
3542-93 
3203-45 
2862-99 
2531-56 
2216-68 
1923-59 
1655-50 
1413-96 
1199-19 
1010-45 
846-298 
704-860 
584-005 
481-532 
395-238 
323-039 
262-975 
213-283 
172-372 
138-849 
111-495 
89-267 
71-270 
56-751 
45-076 
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a {(D 4 * 1 y — i^] A.J. Veyins. 


i 

jp = 0 

p^l 

p = 2 

P = ‘d 

0 

999247 

46-3548 

178-7689 

615-388 

1 

8-87162 

16-3419 

177-6337 

613-923 

2 

7-38673 

44-4156 

175-2606 

608-892 

3 

5-95414 

40-9800 

170-6122 

600-643 

4 

4-704(11 

36-6645 

163-2380 

588-5-20 

f) 

3-66708 

32-0061 

153-3097 

571-636 

G 

2-83015 

27-3902 

141-3932 

549-461 

7 

2-16779 

23-0619 

128-2057 

522-054 

8 

1-65057 

19-1572 

114-4568 

190-0-25 

9 

1-25071 

15-7336 

100-7589 

45t-380 

10 

-943940 

12-7974 

87-5890 

416-333 

11 

•71O(»30 

10-3226 

75-2842 

377-131 

12 

•532556 

8-‘i6G18 

64-0535 

337-940 

19 

•398458 

(vr.7744 

54-0018 

299-757 

U 

•297484 

5-20435 

45-1531 

263-374 

15 

•221679 

4-(l9735 

37-4736 

229-369 

16 

•164913 

3-21139 

30-8897 

198-119 

17 

•122500 

2-50685 

25-3054 

169-824 

18 

•0908724 

1-94975 

20-6139 

144-540 

19 

•0673285 

1-51142 

16-7052 

122-211 

20 

•0498295 1 

1-16808 

13-4733 

102-699 

21 

•0368414 

•900236 

10-8190 

85-8106 

•10 

1 

•0272136 

•692034 

8-65228 

1 71-3175 

23 ' 

•0200849 

•530730 

6-89346 

58-9776 

24 

•0148120 

•406135 

5-47294 

48-5458 

25 

•0109154 

•310161 

4-33096 

39-7849 

26 

•00803856 

•236420 

3-41683 

32-4717 

27 

•00591623 

•179894 

2-68796 

26-4008 

28 

•00435171 

•136657 

2-10891 

21-3874 

29 

•00319917 

•103652 

1-65045 

17-2666 

30 

•00235068 

•078504 

1-28460 

13-8970 
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a' {(D -f 1)^ — -f 1)^} Bp\ % positive. Venus. 


71 

p = 0 

p-l 


0 

447-951 

f)197-2« 

36045-4 

1 

424-582 

5055-43 

35471-2 

2 

387-464 

4804-84 

34419-0 

3 

343-3r)0 

4470-38 

32933-5 

4 

297-290 

4080-18 

31081-8 

5 

25'2 063 

3660-99 

28946-7 

6 

211-488 

3235-50 

26618-6 

7 

174-796 

2821-32 

24186-4 

s 

142-934 

2431-00 

21731-0 

9 

115-818 

2072-54 

19321-6 

10 

93-1103 

1750-22 

17012-9 

11 

74-3419 

1465-46 

14845-2 

12 

58-9993 

1217-63 

12845-7 

13 

46-5732 

1004-69 

11029-4 

14 

36-5893 

823-755 

9401-82 

15 

28-6228 

671-520 

7960-90 

16 

22-3046 

544-536 

6698-96 

17 

17-3203 

439-428 

5604-44 

18 

13-4070 

353-028 

4663-47 

19 

10-3476 

282-448 

3860-91 

20 

7-96501 

225 117 

3181-42 

21 

6-11593 

178-788 

2609-92 

22 

4-68543 

141-525 

2132-23 

23 

3-58195 

111-684 

1735-17 

24 

2-73298 

87-8809 

1406-89 

25 

2-08143 

68-9647 

1136-77 

26 

1-58250 

53-9832 

1 

915-517 

27 

1-20125 

42-1554 

735-058 

28 

-910479 

32-8454 

588-447 

29 

-689121 

25-5372 

469-778 

30 

-520888 

19-8153 

374-056 
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a [{D -f ly ~ + 1)*) i negative. Venm. 


— i 

p = 0 

p^l 

p^2 

0 

447*951 

5197-28 

36045-4 

1 

450-251 

5215-18 

36119-2 

2 

430-208 

5104-60 

35090-2 

3 

395-032 

4877-77 

34770-1 

4 

351-738 

4558-96 

33395-3 

5 

305-713 

4177-06 

31627-6 

6 

260-619 

3760-12 

29548-7 

7 

218-701 

3332-30 

27250-1 

8 

181-139 

2912-63 

24823-7 

9 

148-386 

2514-81 

22355-1 

10 

120-4216 

2147-80 

19917-4 

11 

96-9397 

1816-59 

17569-8 

12 

77-4986 

1523-12 

15356-6 

13 

61-5605 

1267-08 

13308-0 

14 

48-6395 

1046-62 

11441-6 

15 

38-2436 

858 965 

9765-07 

16 

29-9390 1 

700-836 

8277-61 

17 

23-3457 

568-756 

6972-43 

18 

18-1396 

459-303 

5838-58 

19 

14-0489 

369-237 

4862-34 

20 

10-8485 

295-593 

4028-67 

21 1 

8-35441 

235-725 1 

3321 98 

22 

6-41766 

187-307 

2727-03 

23 

4-91853 

148-339 

2229-22 

24 

f 3-76152 

117-111 

1815-13 

25 

2-87098 

92-1892 

1472-47 

26 

2-18721 

72-3729 

1190-34 

27 

1-66343 

56-6708 

959-094 

28 

1-26303 

44-2685 

770-368 

29 

-957552 

34-5027 

616-958 

30 

-724926 

26-8335 

492-720 
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a' {(Z) + !)>“- t-'lC/. Venus. 


i 

p = 0 

p^l 

0 

9185-44 

163970- 

1 

9065-52 

162567- 

2 

8722-88 

158461- 

3 

8199-80 

151939- 

4 

7548-99 

143425- 

5 

6822-81 

133414- 

6 

6067-08 

122420- 

7 

5318-32 

110922- 

8 

4603-33 

99340-4 

9 

3940-00 

88020-2 

10 

3338-69 

77225-0 

11 

2803-95 

67142-4 

12 

2336-00 

57890-6 

13 

1932-08 

49530-5 

14 

1587-52 

42077-2 

If) 

1296-63 

35510-7 

10 

1053-26 

29786-5 

17 

851-310 

24843-7 

18 

684-906 

20611-9 

19 

548-705 

17016-9 

20 

437-852 

13984-5 

21 

348-138 

11443-0 

22 

275-855 

9325-72 

23 

217-905 

7571-41 

24 

171-617 

6125-22 

25 

134-788 

4938-65 

26 

105-589 

3969 32 

27 

82-5145 

3180-71 

28 

64-3347 

2541-55 

29 

50-0521 

2025-38 

30 

38-8608 

1609-93 
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Jupiter. 


i 



a''A^ 

0 

2-018805 

•0385392 

•0209301 

1 

•194930 

•200512 

•0086394 

2 

•0281439 

•0571840 

*0304292 

3 

•00451137 

•0136850 

•0140702 

4 

•000759075 

•00306243 

•00467341 

5 

•000131349 

•000661353 

•00133908 

0 

•0000231 172 

•000139706 


i 




0 

•118302 

•490177 

•116115 

1 

•118939 

•255013 

•163696 

2 

•0284421 

•0891538 

•0991362 

3 

•00030279 

•0262769 

•0419796 

4 

-00137387 

•00704369 

•0147155 

5 

•000290190 

•00177713 

•00159230 

0 

•0000003839 

•000430117 

•00132477 

i 

a"CV 

a"CY 


0 

•0930105 

•228991 


1 

•0418422 

•139507 


2 

•0137512 

•0589828 


3 

•00391940 

•0*206353 


4 

•00102857 

•00642907 


5 

•000255794 

•00185218 


0 

•0000012487 

•000504318 


i 




0 

•0219945 



1 

•0128071 



2 

•00525704 



3 

•00180378 



4 

•000554293 



5 

•000158071 



0 

•0000427049 




[sect. 


a"A.^ 

•00174900 

•00333975 

•00192985 

•00522212 

•003-25231 

•00137257 


•0540104 

•0402*207 

•0473010 

•0325014 

•0100281 

•00048290 

•00230235 
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Mars. 


i 


a"JV 

0 

229114 

4-49988 

1 

•804571 

3-75928 

2 

•405593 

2-89074 

3 

•224607 

2-13691 

4 

•129983 

1-54309 

5 

•0771836 

1-09704 

C) 

•0466132 

•771258 

1 

•0284901 

•537681 

8 

■0175701 

•372392 

9 

•0109110 

•256560 

10 

•00681334 

•175994 

11 

•00427387 

•120293 

12 

•00269104 

•0819091 

13 

•00169982 

•0557079 

14 

•00107664 

•0377737 

15 

•000683546 

•0255G19 

If) 

1 

•000434874 

0172673 

1 

17 

•000277173 

•0116459 

18 

•000176948 

•00784346 

19 

•000113128 

•00527580 

•20 

•0000724205 

•00354460 

21 

1 -0000464159 

•00237898 

22 

•0000297811 

•00159513 

23 

•0000191267 

•00106801 

^^4 

•0000122950 

•000715304 

25 

•00000791006 

•000478450 

26 

•00000509286 

•000319801 

27 

•00000328134 

•000213621 

28 

•00000211557 

•000142609 

29 

•00000136480 

•0000951489 

30 

•000000880963 

•0000634503 


i 




0 

•805995 1-07364 

19'0483 

55*6199 

1 

1-228084 1-02186 

18*4080 

55*2557 

2 

1-050895 1-23434 

16*5922 

53*5762 

3 

•814457 1-30473 

14-2044 

60-3183 

4 

•604237 1-24412 

11-7054 

45-7478 

5 

•437186 1 -10620 

9-37056 

40-3654 

6 

•311182 -936535 

7-33356 

34-6869 

i 



a''Ai 

20 

•00150106 -0148967 

•0943680 

23 

•000453886 00518273 

•0379559 

26 

•000136149 -00175826 

•0146198 

; 




20 

•0827974 

•940514 


23 

•0281560 

•361317 


26 

•00938295 

•134286 


i 

0 

16-3697 



1 

15-5930 



2 

13-HG35 



3 

11-7387 



4 

9-58972 



5 

7-62374 



6 

5-93292 



20 

•0649997 



23 

•0220512 



26 

•00733458 
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Mercury. 



i 



a' A. ^ 


a'A^ 

0 

2081981 

•179759 

•1252855 

•0440444 

•0241364 

1 

•411140 

•464376 

•0914646 

•0528043 

•0227489 

2 

•120179 

•257705 

•167519 

•0476530 

•0255266 

3 

•0389015 

•122616 

•138789 

•0714497 

•0252692 

4 

•0132039 

•0548838 

♦0890266 

•0719791 

•0335508 

5 

•00460G55 

•0237688 

•0502928 

•0563902 

•0371389 

6 

•00163624 

•01008*26 

•0263323 

•0379722 

•03,33614 

7 

•000588587 

•00421654 

•0131072 

•0231604 

•0257195 

8 

•000213722 

■00174517 

•00629291 

•0131806 

•0177725 


1 


rt7)V 


n'B./ 

n'B,^ 

0 

M11682 

1977295 

1-772252 

1-6071 13 

1-259825 

1 

•610090 

1-626068 

1834149 

1-536827 

1-26,5314 

2 

•289401 

1-048.542 

1-589323 

1-501504 

1-214326 

3 

•129374 

•595249 

1-164074 

1 -,3462 14 

1 •16,3209 

4 

•0559922 

•312846 

•756591 

1-072807 

1-062710 

5 

•0237426 

•156179 

•451861 

•770384 

•894799 

6 

•00992677 

•0751578 

•253663 

•508904 

•691321 

7 

•00410781 

•0.351878 

•1.3.5874 

•314718 

•194563 

8 

•00168656 

0161264 

•0701629 

• 1 84GC0 

•331247 


i 


a'C\* 

«'CV 

a'7)„* 

0 

•768937 

3-0 1626 

5-98,351 

•619377 

1 

•586193 

2-64675 

5-68166 

•,5354.36 

2 

•362211 

1-93221 

4-76074 

•384900 

3 

•200324 

1-24869 

3-56256 

•245667 

4 

•10,3463 

•741711 

2-43259 

•144511 

5 

•0510119 

•414433 

1-54498 

•0801264 

6 

•0243209 

•221116 

•9*26391 

•0424868 

7 

•0113048 

•113801 

•530.331 

•0217562 

8 

•00515116 

•0569059 

•292314 

•0108.334 
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FACTORS 
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Values of Mi and 

of their derivatives. Argmnents in the 



Moon portio'ns of JR. 



Terms 

containing iV in the argument 

are to lie 

multiplied by 


For brevity, 

has been printed n Mi 



Argument 

0 


±2V 

i iiV 


Mr 

+ -99276 

+ 1 0005 

+ 1-259 

+ 1-65 


M, 

- -0140 X -2 

- -0835 

- -201 



M, 

0- 

+ -0107 

+ -040 




- 1-3240 

- 1-3478 

-1-701 

- 2-21 


H ~ M. 
an 

+ -0540 X 2 

+ -3245 

+ -852 

+ 2-0 


•‘i"- 

()• 

+ -0509 

± -302 

0- 


i”- 

+ -0842 

+ -1009 

+ -141 




+ *0558 X 2 

+ -2474 

+ -377 




0- 

+ -0368 

± -029 



u 

</k 

- -5390 





±M 

rfk 

+ -0163 X 2 





Argument 




i(/ + 20 

i0-2r) 


- -49080 

- -3452 

- -0736 

- -370 

- -978 

1 i/, 

- -2862 

- -421 

- -799 

- -73 

- 1-10 


+ -3108 

+ -477 

+ *889 

+ -83 

+ 1-34 

n d .. 

- — i/j 

e dn 

+ -5926 

+ -1813 

+ 1-1544 

+ -065 

+ 1-901 

11 d ,, 

- , T/a 
e dn 

+ -7523 

+ 1-1117 

+ 2-5504 

+ 2-151 

+ 3-595 

7L d 

- -j- 

e dn 

+ -8423 

± 1-3085 

+ 2-9124 

+ 2-509 

+ 4-643 



Mi with a sufficient approximation 
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AUXILIARY NUMERICAL TABLES 


[sect. 


Argument 


±(2«+r) 

± (9d - V) 

- *0618 

- 0247 

- -1078 

Im, 

- -2069 

- *24 

- *73 


+ -2065 

1+ 

+ *73 

dn ^ 

4. -0986 

- -001 

+ -230 

^ dn 

+ ‘Mb'd 

+ -580 

+ 2-(> 1 4 

n d 
dn ^ - 

± *5441 

+ -582 

±24)01 


e j — with a suthcient appioxiiiiatiou 


Argument 

M, 

M, 

M, 

±21) 

- *00701 
+ *9870 
± *9868 

± (21) + V) 

- *00291 

- 1*1735 
+ M739 

±(2D-0 
- *03693 

+ 3*1571 
±3*1561 

±(2D + 2Z') 
- *0042 
+ *0826 
± -0818 

±(2l) -2Z') 

- *1265 
+ 7-1584 
±7*1534 

dn 

+ *0254 

+ *0103 

+ *1391 

+ *017 

+ *500 

d 

dn 

- 1*3144 

+ 1*6697 

-4*3211 

- *239 

-9*928 

dn ^ '* 

+1 

± 1*6709 

+ 4*3169 

1+ 

+ 9*892 


- *0052 

- *0013 

- *0214 

+ *009 

- *061 


~ *1445 

+ *625 

- *4104 

+ *004 

- *838 

d.^-' 

+1 

! - 

1 CU 

± *625 

+ *4100 

± *004 

+ *836 

Argument 

±(/-2D) 

± (f - 2D + 1') 

±{i-2D-r) 

± (Z - 2D + 21') 

±(Z-2D-2Z') 


- -08582 

- *2948 

- *0520 

- *7197 

- *0716 

i-"' 

- 1*4897 

f 

- 4*321 

+ 1*188 

-9*11 

- *03 


± 1*4944 

± 4-349 

+ l*i78 

±9*21 

± *05 


4 *2220 

+ *7383 

+ *2751 

+ 1*785 

- *425 

n d 

— y A/g 

e dn 

+ 1*8274 

+ 4*797 

- *501 

+ 9*01 

+ *24 

” ^ if 
e dn = 

+ 1*8518 

+ 4*941 

± *440 

+ 9*67 

+ -31 




= M, 
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VALUES OF MOON FACTORS 


Argument 

Im. 

±(2i-2D) 



± {21-2D + V) 

± (21-2D-1') 

+ *0758 



+ -2552 

+ -0785 


+ *6487 



+ 1-73 

- -18 


+ -5969 



+ 1-53 

± -32 

''' \i 

da ‘ 

- -2228 



- -663 

- -400 

n d 
f" da ^ 

- -8004 



- 1-81 

- -93 

M, 

dti 

± -6010 

d 

did 


± -93 

± -27 


e 

Ml 

- 2 J/* 


Argument 

1 2F 



i{2F + 0 

d. (2F - 1') 


+ 2-9901 



+ 2-8805 

+ 3-0882 


- -2001 



- -25 

- -37 

i-'/3 

+ -1717 



+ -19 

+ -27 

“ J/ 

k'^ dn ‘ 

- 2-9549 



~ 2-937 

- 2-939 

k '' dn 

+ -3785 



+ -509 

+ -645 

n d .. 

9 d u 

± -2873 



± -275 

± -349 


k 

d 

d\i 

Mi 

- 2Mi 


Argument 

± (2F - 2D) 



± (2F - 2D + V) 

-t.(2F-2D-r) 


- -2062 



- -8051 

- -1653 

i 

+ 1-9756 



+ 6-22 

-2-31 


+ 1-9690 



+ 6-18 

±2-31 

- j/ 

k'-' (/« ‘ 

+ -4637 



+ 1-967 

- -319 

- -- J/., 
c/n 

- 1-9243 



-6-282 

+ 2-583 

n d 

k* dn ' ’ 

±1-9019 



±6-164 

+ 2-605 


d 
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Argument 

\m. 

22 + 2^ -[• h — h!^ 

22 + 2ff + A — h'* + V 

2l + 2o + h-h"~V 

+ 1-0000 

+ -8917 

+ 1-1175 

- 1-.3401 

^ nr 

k dn 

- 1-1646 

- 1-0207 

Argument 

Im, 

h - h" 

h - h" + /' 

h - h" - V 

-1-0004 

- 1-0838 

- -9285 

k dn 

+ 1-16399 

+ 1-3385 

+ 1-0256 

Argument 

2T--h-h" 

2T~h-h" + l' 

27 ’ - h~ h" - V 

C-. 

- -04.57 

- -1887 

- -0041 


+ -1213 

+ -7669 

)) 

- -0592 

Argument 

2D + /i - h" 

2D + ;<- A" + 2' 

2D + h - h" ~ V 

Im, 

4- -0355 

+ -1421 

- -0078 

- u 

k dn ‘ 

- -0857 

- -3682 

j) 

+ -0227 


Argument 

D 

4D - 32 

4D - 4F / - 2F 

32 -2D 

M, 

+ -11321 

-0022581 

•00 1 -3-497 1 

+ -0-20801 

il/.3+ M, 

•t--7182[ai - 

•04854 [ 

- -0.334 [k^ - -.3876 

ek'^ + -04856 [ < 

M,~M, 

- -oiso) 

•00066 ) 

•0004) + •9779) 

- -07420; 

Argument 

4D - 2 - 2F 

D-22 

3D - 2F 

2-D 

M, 

+ -011081 

-•1387] 

+ -09511 

+ 003961 


+ -02282 [ ek‘ 

+ •048 

■e'“a, +1-960 [ k-aj 

- -04627 [ eaj 


-•00056) 

+ •310 J 

- 006 ) 

- -58098) 

Argument 

2 + 2f/ + 2i - h” 

2 + 21 - 

1 

.JS 

+ 

i 

3</4-4A-;i"-27' 


- l-49ek 

+ -lOek — -OHOk* 

- •0154e"k 


Argument 

D 

D-2' 

D + 2' 

2-D 

2 + 2'-D 

if« + M, 

+ 1-972 

+ 5143 ' 

+ •815 1 

- -4620] 

-1-38 1 

M,-M, 

- -0280 

- -134 

- -056 

, - 2-465 

1 e 

-5-68 

y <3 ) 

il/g + 

- -0139 

- -058 

- -389 1 

+ -00511 

1 + -023 


- -000 

+ -0009) 

' + -0005 J 

+ -0227; 

) + -089; 


Argument 3D - 2F 3D - 2F - V 

- *326^ + -62^ 

+ •012[ -0 I 

+ +l-959p -3-48p^ 

-000 J *0 J 
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vj 


9 ; 


Extract from the computations for ilhistration of the method. 

,/, = (/) + (29r-26F-3tjr') 

Leverrier (270)‘-(ll*2), i = 26; K^ = A'‘ 

Lev. coef. - [J (27 + 65i + 42t^ + Si^) + |(17 + 1 7i + 41^) + (5 + 2i) 4/ + A^^ 

/^=:28453yl„+ 1581-5^ + 57 ^.,+ ^3 for i= 26 
JJP =. 30034 - 5^1 + 3277/4 + 174.4 + 4.4 + P^p\ 

From the tables for Venus /*= + 6-568, 7 JP-- + 178-87, {{D + 1)^ P - 82 9-8 - 

Factor Y P 3 = 26 (Z) + 2) P - + 4992. l\ - i /\ + (i) + 2) P + (26'^ - 1) 7^ - + 5041 


f 

e 


log 

-3-1104 

- -01)63 
G-6829 
1-0396 


-27-85 

1 

s 

f 


prod. -2-8392 
gy 7-7707 
fac. -4-6099 


log 

- 2-5552 
2-0068 

- 1-6016 
insen. 


^ insensible 

^ 8 3 

/>, +2-9190 h{P., + P,) +3-7003 hJP^- P,) +1*3892 
J/i - 1-6909 +2-3909 M.i+ M, - 1-7760 

fac. -4-6099 fac. - 4-6099 fac. -4-6099 


+ 1-2198 
+ 0"-166 


-2-7011 

-0"-050 


+ 3-7751 
+ 0"-006 


Portion in this term depending on 

All found from the tables except 
[(i; + l)*-26*]7?3% 
whi«h is only needed within 50°/^ of 
value. It can be computed from the 
formulae, but the following estimate 
with sufficient accuracy : 


M ^’1 


1-8724 
62-82 
E, 1047-9 
J^3(e8t.) 12000- 


its true 
general 
gives it 

A 

■ 53270 

- 99360 

- 59730 
. 12000 


-224000 


etc. 


Replace Ajj' by - \ 

The three parts are 

PxM,, ^{Pr-P,)(M^^M,) 


P 

M 

fac. 

"2 

7 


-5-3502 

-5-7466 

- 4-0055 

- 1-6909 

+ 2-3909 

-1-7760 

-4-6099 

-4-6099 

-4-6099 

4-9429 

"4 9429 

4-9429 

- 2-5939 

+ 3-6903 

-3-3343 

-0"-039 

+ 0"-005 

-0"-002 


BR. 


5 



AUXILIARY NUMERICAL TABLES 


<^ = (^) 4- (29^’ - 26 r- 3w'), factor e | 
Portion in this term with factor V 


Levr. coef. = (272)* - (114), i = 26, K^=:A^ 

( 1 1 4) =■ .^V ( 243 - 858i - 1 3 1 7i^ - 926i» - 288i^ - 32?) 

+ -r}i (93 - 4l6i - 74^*2 - 344i" - 48?) + J- (309 + 163i - 6? - 8?) A,^ 

4 \ (123-4- oU 4- 4^•*) ^ 3 * f (29 4- ei) A/ + 5A,^ 


i-26 

i\ 


-008038 

- 177200 

Resultiii) 

yl; 4 - -2364 

- 280300 

From 1\ — 0"’ 

^2*4- 3'417 

- 79800 


^ 3 ' 4 - 32-47 

4 - 67500 

» 2 ^ ■ 

^ 4 * 4 - 270- (est.) 

4 - 50000 

P 2 -P, 

4- 2000- (est.) 

4 - lOOOO 

” 2 ’ 


-410000 

)‘- 


Portion with factor e j . 

//\r> 

Lovr. coef. is approx. 4 times that for U \ , and J- approx.; therefore take ij of 


We obtain 


from I\ 

” 2 


-0"-004 
+ 0"-001 


<I>-{1 + 1') + {2ST - 26 V - 2ot ), factor e^' . , 

Same method gives from 

I\ 4 - 


, ^term ee' — y"'-* est. - J of principal part,^ 


4-0" 008 
- 0''007 

•000 


"'ftTl'.r") 


Other terms insensible. 

Summary for argument (^) -t- (297'- 26 P- 3tz;') 

P 4- P P P 

^From 

4-0"-166 -0"-050 4-0"'006 

- 39 4- r> - 2 

- 64-2 
2 ) - 4 I 

1, 7 '®) 4-7-7 


4- 0"*124 


- 0"-049 


4-0"-004 = 4-0"-079 



TERMS RETAINED BY THE SIEVE 


The terms and coefficients retained and calculated hy the sieve. 


Argument 

h! 

log 1 C 1 

T 

3548 

-1-2 

47’- 3r 

3110 

-1-4 


2650 

- 1*5 

T-V 

2220 

+ -1 

6T-4:V 

1782 

-1-3 

2T- V 

1328 

-10 

3T-2V 

891 

-1-8 

HT-rjV 

453 

-1-5 

57'- 3 K 

438 

- -7 

137’-8r 

15 

-2-0 

•2/^_87’4-5F 

71 

- -3 

„ -57’+ 37 

820 

-1-7 

•2g + 2h 107'+ eV 

74 

-M 

„ - 57’+ 37 

364 

-10 

(j + h-6T + SV 

37 

- -3 

h-6T+4V 

1592 

- 1-8 

A-3r+27 

700 

-1*8' 

A-82’ + 57 

263 

-1-0 

/, + 57^>.3 7 

247 

- 1-6 

h- 5r+3V 

629 

- -9 

1+ 37’-107 

1-85 

0-0 

1+ IGT-IHF 

13 

+ 1-7 

U2l(T- 7) 

425 

- -3 

Z + 247’-237 

467 

- -9 

^ + 297’- 26 7 

28 

+ -1 

1+22{T- V) 

1795 

- 1-5 

(2D-0+ 8r-127 

87 

- 1*2 

„ +217’-207 

102 

+ 10 

„ +137’- 157 

351 

- -7 

„ + 267’- 23 7 

336 

-1-3 

„ +167’ -17 7 

541 

-12 

„ +187’-187 

789 

- -1 

+247’-237 

993 

-1-7 

„ +237’- 21 7 

121^6 

-1-7 

„ +197’-197 

1431 

-10 
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= (0 + (292'- 26 r- 3®'), factor e 


Portion in this term with factor j . 

Levr. coef. =(272y - (114), i = 26, K^ = A^ 

(1 14) - (243 - 8f)Ht - 1317t' - 926? - 288? - 32?) 

+ (93 - 4 1 6i - 741? - 344i» - 48?) + \ (309 + 1 63i - 6? - 8?) 

-f I ( 1 23 ‘+ 5U 4- 4?) Ai f (29 + 6i) Ai + 


^---26 1 \ 
[(2?4l)‘‘*-?]il,;+ -008038 - 177200 

„ yl/+ -2364 -280300 

’^-417 - 79800 

„ A.^+ 32-47 + 67500 

„ J;+ 270- (est.) + 50000 

„ ilg* + 2000- (est.) + 10000 

-410000 a 


Resulting coefficients. 

From Pi — 0''-006 (- of prin, pt.) 

„ +0"-002 


Portion with factor ® j • \^j • 

/e'\' e"'^ 

Levr. coef. is approx. 4 times that for ( 9 J , and J approx. ; therefore take j of 


coef. for 


We obtain 


from Pi 

P,-^P, 
” 2 


<l) = (l + l') + (2ST — 2 (jV-2w'), factor e<»/ . , pterin est. - J of principal part.^ 


Same method gives 


from Pi 


Other terms insensible. 

Summary for argument (?) + (297’— 26 F- 3®') 


■«' 57 ( 1 , 7 "*) 


+ 0"-l24 


+ 0"-004 = + 0"-079 



TERMS RETAINED BY THE SIEVE 


The terms and coefficients retained and calculated hy the 


Argument 


T 

4r--3F 
4r-2F 
T~ y 
^T-iV 
2T- V 

yr-^v 

8^r-5F 

57^-3r 

13r-8F 

2/, -87^+5 F 
„ -r)r4-3F 
2^ + 2/t -lOr+GF 
„ - 5r+3F 

<7 + 4-5r + 3F 
h-iST+4:V 
A-3r+2F 
A-82^+5F 
;, + r)r~3F 
h-bT^ZV 

1+ yr-iov 

1+ lGr-18F 
/+ 21 (T- V) 
l+2Vr-23V 
U29T-2^V 
1+22{T- V) 
(2D-0+ 8r-12F 
„ +217^-20F 

„ +13r-15F 

„ +267’- 23 F 

„ +167’-17F 

„ +187’-18F 

„ +247’-23F 

„ + 237’- 21 F 

„ +197’-19F 


h! 

1 log 1 C 1 

3548 

- 1-2 

3110 

-1*4 

2656 

- 1-5 

2220 

+ -1 

1782 

-1-3 

1328 

-10 

891 

-1-8 

453 

-1-5 

438 

- *7 

15 

-2-0 

71 

- -3 

820 

- 1-7 

74 

-M 

364 

-10 

37 

- *3 

1592 

-1-8 

700 

-1-8 ‘ 

263 

-1-0 

247 

-1-6 

629 

- -9 

1-85 

00 

13 

+ 1*7 

425 

- -3 

467 

- *9 

28 

+ -1 

1795 

- 1*5 

87 

- 1-2 

102 

+ 10 

351 

- *7 

336 

-1-3 

541 

-1-2 

789 

- -1 

993 

-1-7 

1226 

-1-7 

1431 

-10 


sieve. 






TERMS RETAINED BY THE SIEVE 


Argument 

l«! 

log|C?| 

3D-F + ^+27'-3<2 

100 

- -9 

2F-/ + 3r-4g 

68 

- M 

T~J 

3200 

- -7 

2{T~J) 

6400 

-1-6 

J 

299 

- -5 

T %T 

2900 

^ 1-8 

%T 

598 

- 1-8 

h + J 

108 

-1-3 

h - %T 

789 

- 1-8 

2h-2J 

980 

-1*4 

'2h^J 

681 

1-3 

(2/» + 2j7)-2./ 

204 

- -1 

„ -3J 

93 

- *7 


503 

- -9 

„ -1./ 

395 

- 1-7 

g ^h- J 

102 

- -7 

g^2J 

6-57 

- 1-6 

T~M 

1662 

-1-6 

2T-2M 

3323 

-1-9 

2T-'nT 

1437 

1*8 

T-2M 

225 

- M 

2T-^M 

450 

- 1'8 

h-T+2M 

34 

-1-0 

2g + 2h + yr-^M 

127 

- 1-7 

44 + 57^-94/ 

•80 

- -6 

2g^h + lT~\M 

5-46 

1-4 

/- 267' 4- 244/ 

57 

1 -3 

2D -^-207’+ 164/ 

40 

- *2 

D- 237'+ 204/ 

8-65 

•5 

4D-^-2F-20(7'-4/) 

28 

-1-7 

3D- F-/-207'+ 183/ 

0*04 

- •4> 

4D-3/-237'+254/ 

•58 

+ 1-4 

S 

120 

computed 




SECTION VI. 


INEQUALITIES IN LONGITUDE, LATITUDE AND PARALLAX DUE TO 
THE DIRECT ACTION OF THE PLANETS. 

In all cases <f) denotes the argument of the primary term. 

When a primary term consists of two or more terms of the same period 
but with different arguments, the latter are denoted by <^o + ^ where G 
depends on fi\ ot', and on the coefficients. Each term is set down. The 
single term, argument <j>, which is obtained by combining them, is then 
given. The secondary terms, arguments where <!>' consists of the lunar 

arguments only, follow. Secondary terms having the same arguments as 
certain primary terms have not been added to the latter but are given 
separately in their proper places. (See Addendinn.) 

Coefficients of Sines in Longitude. 


Venus. Short Penod Primaries. 
c}>^i(T- V) 
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VENUS TERMS IN LONGITUDE 
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4, = l+i{T- V) 


+ + i{T-V) 


i 

0 

0 

-/ 

0 

-/4-2i> 

-14 



+ 0' 

•001 



- 13 



4- 

2 



- 12 



4- 

3 



- 11 



4- 

3 



- 10 



4- 

4 



- 9 



4- 

5 

- 

o 

o 

- 8 



4- 

6 

- 

1 

- 7 



4- 

8 

- 

2 

- 6 



4- 

10 

- 

2 

- ^ 

+ 0" 

002 

4- 

12 

- 

2 

- 4 


3 

4 

15 

- 

3 

- 3 


3 

4- 

18 

- 

4 

- 2 

+ 

4 

4- 

20 

- 

4 

- 1 

+ 

5 

4- 

21 

- 

4 

0 

4- 

r)t 


* 

- 

4 

1 

+ 

3 

4- 

17 

- 

3 

2 

4- 

3 

4- 

13 

- 

3 

3 

+ 

2 

4- 

11 

- 

2 

4 



4- 

9 

- 

2 

5 



4- 

7 



6 


i 

4- 

6 



7 


1 

4- 

5 



8 


j 

4- 

4 



• 9 


1 

4- 

3 



10 


i 

1 

) 

4- 

2 

1 



* Argument is zero. 

t Argument is I and this coefficient therefore pro- 
duces a slight change in the value of c, which is 
however too small to affect any term sensibly. 


i 

0 

(p — i 

0 


4- 0"-001 

1 


4- 2 

2 

<D 

4- 2 

3 

‘co 

a 

<v 

4 1 

4 

CA 

4- 1 

5 


4- 1 

G 


+ 1 


-I'-l 

(7^- V) 

i 

0 

<p-l 

- - - 

L__ — 

— 

1 


4 0"-001 

0 


4- 1 

1 


i + 2 

2 


4- 3 

3 

o 

4- 3 

4 

2 

’w 

4- 3 

5 

C 

0) 

w 

4- 3 

G 

.2 

4- 2 

7 


4- 2 

8 


4- 2 

9 


4- 2 

10 


4- 1 

11 


4- 1 
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INEQUALITIES DUE TO DIRECT ACTION [SECT. 


</) = 2Z)-Z-r-i(r- V) (l, = 2D + l + i(T-V) 


i 

<t> 

0+/ 1 

_ i 

+ l — 2D 


i 

0 

<f> — 1 

-1 



+ 0" 

•002 



-13 


o 

p 

b 

f 

0 



+ 

4 


' 

-12 


1 

1 

+ 

1 

+ 

7 

- 0" 001 


-11 


— 2 

2 

+ 

1 

+ 

8 

_ 2 


-10 


- 2 

3 

+ 

2 


9 

2 


- 9 


3 

4 

+ 

1 

+ 

10 

— 2 


- 8 


3 

5 

+ 

1 

+ 

9 

— 2 


- 7 


4 

6 

4- 

1 

+ 

9 

_ 2 


- 6 

.'2 

'co 

4 

7 

+ 

1 

'i + 

9 

— 2 


- 5 

QJ 

W 

.s 

5 

8 

+ 

1 

+ 

7 

1 


- 4 


G 

9 



1 + 

0 

1 


- 3 


6 

10 



+ 

C) 

1 


- 2 


G 

11 



: 

4 



- 1 


G 

12 



' + 

3 



0 


5 

13 




3 



1 


2 

14 



+ 

2 



2 


1 

15 



+ 

i) 














<f) = 2D 

-i-r 

ST 

1 

T 

i 

0 

0 4 i 

2 


- 0"*002 

3 


2 

4 


3 

5 


•> 

“ O 

G 

2 

— 3 

7 


3 ; 

8 


3 

9 


_ 2 

10 


_ 2 

11 


_ 2 







76 INEQUALITIES DUE TO DIRECT ACTION [SECT 


</) = 2Z-2i) + /'+l(2’- V) 


— 

i 

0 

(p -1 

<p-I + 2D 

(p — 21 

1 


- 0"-003 



2 


6 



2 


12 

+ 0''-002 


4 

0"-005 

42 

+ 9 

- 0"*003 

5 


+ 

34 

7 

+ 2 

6 



12 

— 2 


7 


+ 7 



8 


+ 

5 



9 


+ 

4 



10 


+ 2 




4 > = - 2 l-‘in + l" + i(T- V) 


i <p 

<p-l 

<p — 1 21) 

3 

+ 0"-002 


^ S, 

+ 4 


5 •'1 

+ 13 

- 0''-003 

fl 



ft a» 

0 25 

10 

+ 2 

7 

4 


8 

— 2 



4,=^2l-2D-l' + i(T- V) ^=^2l-2D-l" + i(T-V) 



<j> = 2l-2D ±l{i + 2) T-iV} 

Soft long poiio<l terms. 


<f> = h ± \(i + 2}T-iV} 
See long period terms. 




vij 

VENUS 

TERMS IN LONGITUDE 
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Venus. 

Long Period Primaries. 


</)o=8!r— sr, s 

= - 452"-9 <l>„=i:iT-SV, s 

= - 14"*85 

period 7*8 years 

period 239 

years 

^0 — — 2h/ 

+ 0"-007 

00 — to' — 4/t" 

+ 0"-0028 

00 — w" — 2//' 

— 2 

00 — ot'' — \h" 

9 

00 — d*!*! 

+ 1 

00 — 3ot' — 2^" 

+ 31 

00 — 2'ZJT^ — 

1 

I 

1 

> 

1 

30 

4>„+ur-,l> 

+ 0"007 

00 ~ w — 'Iw" - 2//' 

+ 10 



00 — 3rEr'' — 2/i'' 

1 



00 - 39° = 0 

+ 0"-003 


<^o = 

I + ST-IOV, 

= + l"-85 




period 1920 

years 




<^o - 

w' — 6A" 

+ 0' 

'•34 

0 + / 

-t 0"-019 


- Oh" 

- 

5 

0 ± (Z - 2i>) 

+ 4 

<^0- 

1 

> 

CO 

+ 

7 

0+ 2i> 

+ 4 

<^o - 

2tzr' — m" — ih" 

- 

5 



^0 

1 

1 

+ 

1 




b'm' - 2h" 

+ 

1 



^0 ~ 

— vj" — 2h'' 

-- 

1 



<^0_+ 

33° = 0 

+ 0" 

•35 




<f>o=-i- ler+isF, «= 

+ 13"-01 





period 273 years 

00 -2/." -13"-99 

0 + / 

-0" 

•766 

0 + (2/)-^-/') 

~0"-008 

00 - 2w' - 66 

cf>~l 

- 

815 

<l> + 1 — 1' 

5 

(f)Q — m' ~ Tjs" + 71 

0 + 2i) - / 


159 


6 

00 — 2'®’^^ ^0 

0 - 2y>> 4- / 

- 

1(58 

0 + 2/^’ 

+ 29 

00 — XET^ 4- — 2}l' — 8 

0 + (2/y 4- 1) 

- 

20 

0 + 3^ 

4 

00 + — ot" — 2h'' - 27 

0 + 27) 

- 

166 

0±(4/>--2/) 

— 2 

00 + 2rtr' - 4//' — 6 

0-2/> 

- 

168 

0 + (4/>-0 

4 

00 + xer' + ra" — 4/i" + 5 

0 + 2^ 

- 

52 

cfi + {2 P + 1) 

4- 5 

00 + 2sr" — Ml' — 1 

00'=<#> 

0-2^ 

0 ± (2/) - ^') 
<j>±D 

0 + (^ + ^ ) 

4* 

4- 

55 

11 

5 

4 

0±(2/<’-/) 

4- 2 
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INEQUALITIES DUE TO DIRECT ACTION 


[sect. 


= / +29^-26 F, 5 = -27"-85 
period 127 years 

+0"-096 

<^o ~ — 24" — 25 

4>o ~ + 79 

<f>^) — — ■ztr" — 71 

00 — — 2tzr" + 22 

00 — 2 

00 — 2®-' + -m" — 24" + 3 

00 + - 2ot" — 24" — 1 

00 — 470-' + ot" — 1 

00+ 112 ° = 0 

0 ± ^ + 0"’006 


0 = 2Z)--Z+18(7’~ F), 6-= 788"-8 
period 3*5 years 

0 + ^ + 0"’003 

0o = Z-2Z)-8T+ 12F, 5 = + 87"-07 

period 41 years 


00 - 44" -0"-032 

00 — 2'rtr' — 24" — 2 

00 — w' — ot" — 24" + 1 

00-123° = 0 + 0"-0^3 


0o = /-2i)~ I6r+17F, .^ = 539"-9 

period 6*6 years 

f 

00 — TO-' + 0"'004 

00 ■" ~ 2 

00-73° :^0 +0"-003 


0 = / 4- 21 (T - F), s = 425" 0 
period 8*4 years 

0 + 0"-030 

0-/ + 0"-005 

(fi + l + 0"'002 

0= / + 22(r- F), 5 = - 1795" 
period 2*0 years 
0 + 0"-002 

00 — I 2t^T — 23 Y f s — — 465"*8 

period 7-6 years 
00 - tir' + 0"-009 
00 — t«r" — 3 

0„- 88^=0 + 0'^006 

^ = 2i> - ^ + 19 {T- V\ s = - 1430"-7 
period 2 5 years 
0 + 0">002 

</>=/-2i)-13f+15F-24", 6-=-351" 
period 10 years 
0 _ - 0"-025 
-0" 004 


<f>„ = '2D-l + ‘2,lT-%0V, s = -101"-9 

period 35 years 

+0"171 

</>o — ot" — 52 

<l>i, + in' — 2h" — 3 

~ 2k" + 1 

00' = <!> +0"-126 

d>+l +0"-010 

<I>-1 +0"-007 



VI] 


VENUS TERMS IN LONGITUDE 
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00 = 2/ ^ 22) + Sr- 6 7, s = 74"06 

period 48 years 

0O-2A" -0"-036 

<f>Q~ 2m' — 56 

00 — ot'-ot" + 35 

fj>Q — 2m" — 5 

00 + 17'^ 30^ = + 0^^-Q65 

4>-2I> ~ To”-002 

<t>-l +86 

(f>+l + 5 

^ + - 16 

0 — 2/ + 6 

</> - 2/ + 22) - 2 

0 = 22^-22) + 3(f- 7), 6 = 819"-5 

period 4 '3 years 

0_ - Q^^-QQ2 

0-22’ +0"-002 


d!) = 2/ - 22) + 7 r - 5 7 + 1 5°, s = 2296" 

period 1*5 years 

<f> insensible 
d >-/ + 0"^qi^ 

0=2/-22)+9r-..77+2O°, 5=-2145" 

period 1*7 years 
</> insensible 

0-7 -0"*003 


0o = 2/^-2Z) + 6r-5 7, a* = -71"*27 
period 50 years 
0o-'nr' +0"’071 

<;^>o - ot" - 21 

<^o - 90^ = 0 +0^’*Q 54 
0±/ +0"*0b3 

0-22^ -0''-003 


0 = :U - 22)+ 24 (r- 7), 6* = 61"-()7 0 = 42) - / - 22"+ 15 (T- 7), .9=-30"*7 

period 58 years period 116 years 

0 +0"*010 0 -0"*002 


0 

S 

period 

coef. primary secondary 

A-42’+27 + r 

- 2848" 

I 4- years 

+ 0"*005 

/^_57’+3K + r 

- 628"*7 

5-2 years 

+ 0"*016 0~22’ +0"*003 

A- 6^+4^ + A" 

+ 1590*8 

2*2 years 

- 0"-0()9 

*A-72’+5r+ A" 

+ 3810* 

0*9 years 

-^"•003 

0o=/+2^+A+192’— 

207, i'=88' 

'*85 0 = 

l+h+ier-lHV+k", s=-203"-8 

period 40 

years 


period 17 years 

00 + -nr' — h" 

-0"*003 


4> + 0"-008 

00 + — A" 

+ 2 



<^, + 166* = (^ 

+ 0"*002 





80 


INEQUALITIES DUE TO DIRECT ACTION 


[sect. 


</>=^-fiy-f24r-23F-r,5=125"*9</>=Z+4^-f3/i+24r-22F-3r,s=-0"‘07* 


period 28 years 

period 50000 years 


<l> + 0''-003 

+ 0"-02 


-D - nT+UV, s = 50"-0 

(i)=D + 2()(r- V). s 

= - 502"-8 

period 71 years 

period 7*1 years 

^0 — — 24" - 0"*008 

<I>Q — w" — 24" - 5 

- 8'2” ■- <#. +0"’013 

<t> + 0"-002 

i) + 25T-23V, 6 = -64'' y 

= 47^-3V, 

s = 37"-03 

period 54 years 

period 94 years 


<#.„ - 2A" + 0"-006 

+ 12 
<!>(, — — uf ' — G 

<^Q - 2'm" 4- 1 

— lir' — 0"'059'j 

“ zo’" + 18 

+ 'uj' — 24" -- 1 j 

c^o + 93° +0"‘044j 

from term 
of Ji fac- 
tored by — , 

</)(( + 190 ~ <f) + 0 '013 

+0"'005i 

</>o — ot" — 0"'002 

1 from first 
part of Jii 

: W - 2 F+ vrr- 18 V, ,5 = ()" ;33 

period 560 years 

<#>„-88° +0"-004, 

+ 0"'003 

<^„ + 93“-</. +0"-040 


<^0 - zir" — 0"'001 

=^<f> +0"'002 

(f> — 1 +0 '029 

<j) -i- 1 +0 ‘002 

+ -0"'005 

<#> - 2? + 0"'002 



^ = h — h" 

Venus -0"-0152 

Mercury - 1 

Jupiter + 24 

Alars + 3 

-0"-013~ 


* This period is approximate only. But the coefficient is insensible to observation whatever 
be the period. 



VI] 


JUPITER TERMS IN LONGITUDE 


81 


Jupiter. Mort Period Primaries. 


i 

<P 


0 

See long period primaries 


1 


-0"-002 

2 

-0"-002 



i 

0 


1 

-0'''069 

-0"*008 

2 

-0"-001 

-0"‘005 


(l> = W + i{T-J) 


i 

0 

1 

-0"'002 

2 

- 22 

3 

- 10 

4 

3 




i 

0 

</) - 1 

-1 


+ 0"-003 

0 


* 

1 


+ 0"-003 

2 

+ 0"*002 

+ 7 

3 

+ 1 

+ 3 


* The argument is zero. 


^ = l-lD + i(T-J) 


j> = 2l-W + i(T-J) 


i 

0 

0-/ 


<l> 21 

0 

-1 

-2 

-3 

-4 

insensible 

See lo: 

insensible 

-0"-003 

7 

ng period pr: 

+ 23 . 

+ 4 

+ 0"’001 

maries 

-0'''005 

+ 0"'002 


i 

0 

(l)-l 

(t)-l+2D 

0 


-0"-002 


-f 


4 


-2 

-0"'007 

- 39 

+ 0"'008 

-3 

3 

- 17 

+ 3 

-4' 


5 



4,=2l-W+l"+i{T-J) 


“:n 

t 

0 


-3 

insensible 

See long peri 

! 

+ 0"-006 
iod primaries 


BK. 


<f,= W + l+i(J-T) 


i 

0 


-2 

-3 

-0"‘001 

-r-oos 

-0"*002 


6 
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INEQUALITIES DUE TO DIRECT ACTION 


[sect. 


Jupiter. Long Period Primaries. 


(^0 = /, s=29r*13 

4>q — 2f/ , s 

= 598"*3 

period 11 *8 

years 

period 5-9 years 

4)0 — VT 

- 0"-208 

<j>0 - 2ot" 

-0"*008 


+ 0"*021 

(fyQ — vt' — w" 

4- 0"-001 

^,,4- 173 50' = 0 

+ 0"-209 

4,„- 2h" 

-0"-001 


4 - 0"-021 

<#>„ + 162° = 4> 

4- 0"-009 

4> i 

+ 0"*002 



2l-2D-2{T-J),s = - 203''o8 

<f>„=2l-2D+3J 

-2T.s=d5"- 

period 15 ‘5 years 

period 37 years 

4, -0' 

'•187 

<Ao - 

-0"-190 

(j) — 1 + 0 

'•811 

fjio — ■n: 

4-0"-017 

(f> + 1 — 

7 

<^0+173'’ 16' --=<^ 

4-0"*190 

^ 4- 2D — 

•170 

<j> — 1 

4 - 0"'306 

^-.1-.2D + 

7 

4> + l 

4- 5 

<f} — 21 + 

56 

</>-/5 + 2Z> 

58 

4, -21 + -ID 

15 

4,-l-2D 

4- 3 

^ — 1 — 1! — 

4 

4>-2l 

4- 20 

4^-.1-V + 2D - 

8 

4>->2D 

4- 6 

4> ^ 1 + 1' + 

5 

<I> + 2D 

- 2 

4» — + 

4 

4>-2l + 2D 

6 

4>+2D 

13 

4>-l-l'^2D 

3 

<ji — 2D + 

9 

4>—l + l' 

4 2 

<l>-2l + iD - 

2 

4,-l-l'+-2D 

3 



4> — 1 1' 

4 - 2 

21 - 2D + J-2T+^'‘, s=-502"-6 

4>=2l-2D+4J-2T- 

-2^^ 6* =394^ 

period 7 years 

period 9 years 

4> 

insensible 

4, 

-0''-002 


~0"007 

<f> — 1 

-0"-009 


• 

4>-l + 2D +0"-002 

2D-2F+2(T- 

A s = -979"-8 

<j>=h+j-h'', 

s=108"-36 


period 3* 6 years 
< f > + 0^'‘002 


period 33 years 
<#) -0"-004 
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MARS TERMS IN LONGITUDE 
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Mars. Short Period Primaries. 


i 

<t> 


1 

1 

o 

o 

-0"-002 

2 

8 


3 

3 


4 

2 



i 

<^o- w" 

00 -ot' 

-1 

-0"'002 


0 

5 


1 

55 

+ 0"-008 See long period 

2 

+ 9 

— 1 primaries 

3 

+ 4 



= 2i/ 4- ^ (/> = 3# + 3 (If - T) ~ 3t!r' 


i 

00 — 2 “nr" 

00 - tor' — tp-" 

0-0o + 6H« 

1 

-0"’002 


1 

o 

o 

o 

u> 

2 

10 

+ 0"*003 

12 

3 

+ 4 

1 

+ 5 

4 

+ 2 


+ 2 

5 

+ 1 


+ 1 


i 

0 

3 

4 

-0"-002 

+ 1 




i 

0 

(f>- 1 

(p — l-\- 2D 

2 


-0"*002 


3 


5 


4 

-0"-002 

+ 17 

CO 

o 

o 

o 

I 

5 


+ 3 


6 


+ 1 



</>o=2^-2i) + i(T-#) + if 


i 

00 ~ ^ 

00- w' 

0O + 21P-? 

-4 

-6 

-6 

-0"*014 

+ 0"002 

-0"-003 

See long period primaries 

+ 0"-003 


6-~.2 
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INEQUALITIES DUE TO DIRECT ACTION 


[sect. 


Mars. Long Period Primaries. 

<l>o = 2M-T,s = 224"’8 <^o = 2^ - 2Z) + 6ilf - ST, 5 = - 127"*31 . 

period 16 years period 28 years 


4>o 

~0"-055 

<#>o “ 

-0"‘014 


+ 8 


+ 2 

.^„+213' = (^ 

+ 0"-060 

<^o + 211° = 

+ 0"015 


+ 0"*007 

€f> — l 

-0"'038 



<l>-‘l + 2D 

+ 8 




3 

<f>. 

,= 2l-2D + SM - 

er, s = + 97"-56 



period 38 years 



<t>o — 2ts-" 

-0"-015 



— W — TU 

+ 5 




1 



</>„ + 63° = <#> 

-0"‘019 




-.0"-030 



<I>-U2D 

+ 6 



</.-2/ 

- 2 



<f> + 1 

1 


4i)-3/ + 25iti- 

2ST, s = -0"-58 

^0 ^ h 2ili - 

T, s = 34"-07 

period 6000 years 

period 104 years 

<^o - 

+ 0"-03 

^0 — 'ZET — h 

-0"-015 

— vj' — W*' 

1 

"" 737 ■" h 

+ 2 

113“ = <!> 

-0"-04 

<^,,+ 165° = <l> 

+ 0"*017 

<i>±i 

-0"-002 




Mei'cury, No Short Period Pnmaries. 

Long Period Primaries. 

<f>, = l-2D-T + SQ,s = - 90"*36 6T -^Q, s - 449"‘3 


period 39 

years 

period 7*9 years 

<^o - 

-0"*047 

— — 0"*004 

<^0 nr — nr 

+ 3 

</>o - nr' + 1 



40 

(#>, + il3°=<^ +0"-003 

<^o+ 75° = 

+ 0"*075 


4 > — 1 

+ 0"-006 




vi] VENUS TERMS IN LATITUDE 85 

<l>, = -2F+l-8T+4,Q,s = 67"-76 

4> = l + 2g->rZh-T- 3Q+r, s=-100''-4!l 

period 52 years 

period 35 years 

- 0"*004 

<#. + 0"-002 

fl>Q —W 4- 1 


<^„ + 113° = <#> +0"-003 

Saturn. 

Short Period Prim allies. 

Long Period Primaries. 

4> = 2D-1-2{S-T) 

<^„ = S. s=120" 43 

<f> insensible 

period 30 years 

<I> + 1 +0''*003 

^0 — 0"‘026 

(/> = 2Z)-2(^-T) 

<#>o — + 2 

</.„ + 90“ = + 0''-024 

<f> -0"-001 

<f> + l + 0^ *003 

cl)^2l-2D-2{T-Sl s=-5G0"*9 

<l>^2l-2D+2S--2T-w'\ 5=-440"*5 

period 6*3 years 

period 8 years 

insensible 

insensible 

+0"-014 

+0"*004 

<h-l + 2D - 3 




Coefficients of 

Sines in Latitude. 



Short Period Primaries. 



V) 

(/> = 2i) 


F) 

i 

<f)it P <^ + 2^+/ 

i 

^ -^F +1 


1 

+ 0"*005 + 0"*002 

-4 

+ 0"*002 


2 

+ 3 

-3 

+ 3 




-2 

4- 3 


(/> = 2i^ 

-2i)-b^(T- V) 

-1 

+ 3 


i 

<h-F 

0 

4- 2 


0 

- 0'' 003 




1 

G 

<^ = 2Z- 

2D + i{T- 


2 

11 

% 

<t>±P -1 


3 

See long period primaries 

1 

- 0"*002 


4 

+ 25 

2 

6 


5 

+ 9 

3 

See long period primaries 

6 

4- 5 

4 

4- 4 


7 

+ 3 

6 

4- 2 


8 

4- 2 






86 INEQUALITIES DUE TO DIRECT ACTION [SECT. 


=^h + h" 

-iT+(i-2)F 

</> = 2F- 

■W + T+iiT-V)-90° 

i 

<^ + F 

i 

<f> — F 

-1 

+ 0"'003 

4 

+ 0"-002 

0 

+ ’ 4 

5 

See long period primaries 

1 

+ 6 

6 

-0"-002 

2 

+ 8 



3 

+ 10 



4 

+ 18 

(f) =z h ~ 

h" + i{T^ V) 

5 

See long period primaries 

i 

<I> + F 

6 

25 

-7 

+ 0"-002 

7 

9 

-6 

+ 3 

8 

5 

-5 

+ 5 

9 

3 

-4 

+ 6 

10 

2 

-3 

4- 9 



-2 

4- 14 



-1 

+ 27 



0 

See long period primaries 



1 

15 


2-6 

3-3 


iMg Period Primaries. 


</) = 2/ - 27) + 3 (T- F), 5=-3():r-9 
period 9^ years 

-0"-003 

fj) ^ P — I + 32 

+ F — 1 + 2D - 6 

tli-F-l + 2D - 5 

</)-F-2/ + 4 

</)=2F-2D+6r-5F-90°, s=-7r*: 
period 50 years 
<I,-F +0"-068 

<I>~F+2D - 2 

<ji — F — I + 4 

^-F-^l - f 4 

<;> = /i-/i",s = -190"*8 
period 18-6 years 
<^ + /’ -0"-241 

tj) + F — 2D + 8 

<ji + F I — 13 


^ = 2/ - 2 77 + 8 r- 6 F+ 1 7 " 30', s= 74"‘06 

period 18 years 

^ + 7^ + 0^^'003 

<l> + F — I + 4 

<^ = 2F-27) + 3(T-F), ,9 = 819"-5 

r period 4^ years 

<f>-F -0"-044 

<t,~F-l - 2 

<l>-F+l 4 - 2 

- 5r + 3 F, s = - 628" *72 

period 5 6 years 

<I> + F +0"*072 

4> + F-2D - 2 

<^ + 7?’ + ^ + 4 

«/» + 7^-^ - 4 



VI] JUPITER TERMS IN LATITUDE 87 

<^ = Z + 10F+ 33°, + r-85 (^=-Z-16r+18r-150°69', s = 13''*01 



period 1920 

years 

period 273 

years 




+ 0"-016 

<f> ^ P 

-0" 

•650 




<l> ± (P- 2D) 


22 


= 1 + 29r- 267+112'’, s = - 27"-85 

± (F+ 2D) 

- 

12 


period 127 


<l>+ P+l 


79 


years 

F^l 

•_ 

56 


P 

+ 0"'00r) 

<f> + P—l 

- 

7 




th-P+l 

+ 

14 

<#>= 

= 2Z)-Z+2iy-20F~ 

.87°0>‘=-10r-92 

tf>+P-l + 2D 

- 

13 




<j> — P + 1 — 2D 

- 

17 


period 35 

years 

^ + P ■¥ 21 

- 

8 


<f>+ P 

+ 0"-006 

<l> + P-2l 

- 

6 



<l,±{P^l-^2D) 

- 

3 


Jupiter. No Short Period Primaries. 

Long Period Primaries. 

5 = 299"! cl>^2l^2D-2{T^J),s^- 203"-58 

period 12 years period 15*5 years 


<f> + P 

+ 0"‘002 

4>±P -0"*008 

<l,=2l-2D+3J-2T+l7rm',s=d5"-a6 

+ 35 

<f> + P-~l + 36 

period 37 

^ P 

<f> — 1 + P 
cf> - l-P + 2D 

years 

+ 0"009 
+ 14 

2 

<l>--P-l+2D - 6 

<f> + P -1 + 2D - 7 

^-P-2l 4- 4 

^JrP-2l 4 - 2 

<l> — l + P+ 2D 

3 • 


^:=h-h"-]-J-vT",s= 108"-36 

<j[, = /i~/i'^s=-190"-8 

period 33 

years 

period 19 years 

<l> + P 

- 0"005 

^^P + 0"*038 

<I>=^2F--2D-2(T 

- J),s = 979'^-8 

^ — h — 1i T — J^ s = 3058^ 

period 3*6 

years 

period 1*2 years 

tj) — P 

- 0"*026 

4> + P - 0"-002 

<;>-2F~2i)-3(T- 

. J), s=-2269" 


period 1*6 

years 


<h-P 

+ 0"'005 
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INEQUALITIES DUE TO DIRECT ACTION 


[sect. VI 


Mars, 


</) = A + 2A/ - T + 165°, ,9 = 34"*07 
period 104 years 
4> + F - 0"010 


s = -190"-8 
period 19 years 

+ 0"004 


Mercury and Saturn, None. 


Coefficients of Cosines in Parallax. 


0 = 2/-2i) + 3(r- V) 

<f> — 1 

+ 0"'006 

^^/ + 167^-18F+15r 


+ 0"-007 

<f> = 2l-^2D-2{T~J) 


+ 0"-007 

<t>^2l-2D+3J-2r+ 173° 

<fi — 1 

+ 0"*003 


Additions to Annual Mean Motions. 

Perigee + 2''*69 

Node - 1"-42 



ADDENDUM. FINAL VALUES. 


The primary coefficients in longitude having short periods can be combined 
with the secondary coefficients having the salme arguments in many cases. These 
combinations have been made in the following tables which therefore represent 
the final coefficients due to the direct action of the planets, with the following 
exceptions : — 

Venus. All long period primaries together with their secondaries except those 
arising from the arguments 

2l--2I) + 2r + i(T~V) + a, 2F-^2D+3(r- F), 2F~ 2D + 67^- 5F- 90% 

which have been included. 

Jfars. The long period primaries with arguments 

4D-3U25J/- 23^-^113% h + 2M^T+16r)\ 

Saturn and Mercury. All terms. 

These terms, omitted here, can be taken immediately from the tables in Sect. vi. 


Final coefficients of sines of i{T- V) arguments at top of each column. 





2r+26°-6 


2D 


2D- 

-T + 

78“ 

i 


0 

- 4 

4 0' 

'•001 

-15 

+ 0" 

O 

o 

- 15 

4 - 0" 

•001 

1 

+ 0" 

•480 

-3 

+ 

2 

-14 

4 - 

2 

-14 

4 - 

1 

2 

+ 

•200 

-2 

+ 

2 

-13 

4 - 

2 

-13 

4 - 

2 

3 

4 - 

92 

-1 

+ 

3 

-12 

4 - 

2 

-12 

4 - 

2 

4 

4 - 

60 

0 

+ 

4 

-11 

4 - 

3 

- 11 

4 

3 

5 

+ 

38 

1 

+ 

6 

- 10 

+ 

5 

-10 

4 - 

4 

6 

+ 

2.5 

2 

+ 

8 

- 9 

4 - 

6 

- 9 

4 - 

4 

7 

+ 

17 

3 

+ 

37 

- 8 

4 - 

8 

- 8 

4- 

4 

8 

+ 

12 

4 

- 

8 

- 7 

4 - 

8 

- 7 

4 - 

5 

9 

+ 

8 

5 


3 

- 6 

4 - 

11 

- 6 

4 - 

5 

10 

+ 

6 

6 

- 

4 

- 5 

4 - 

11 

- 5 

4 - 

7 

11 

+ 

4 

7 

- 

1 

- 4 

4 - 

10 

- 4 

4 - 

7 

12 

+ 

1 




- 3 


36 

- 3 

4 - 

4 







- 2 

4 - 

26 

- 2 

4 - 

4 

T 

-87' 

’•8 




- 1 

4 - 

15 

- 1 

4 - 

3 

~3 

-0" 

o 

o 




0 

4 - 

16 

0 

4 - 

2 

-2 

- 

4 




1 

4 - 

15 

1 

4 - 

1 

-1 

- 

8 




2 

4 - 

8 




0 

_ 

15 




3 

4 - 

4 

2D- 

2T- 

18° 

1 

- 

47 




4 

4 - 

4 

-6 

-0" 

•006 

2 

+ 

76 




5 

4 * 

4 




3 

+ 

21 




6 

+ 

3 




4 

4 

12 




7 

4 - 

3 




5 

+ 

7 




8 

4 - 

3 




6 

+ 

6 




9 

4 - 

2 




7 

4 - 

4 




10 

4 - 

1 





8 + 1 
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ADDENDUM. 

FINAL 

VALUES 



1 

4- 27^4- 27" 

2D-. 

2-^4- 89" 


21 

-8 -0"*002 

2 

4- 0"-001 

-10 

- 0"*001 

-3 

- O'' -001 

-7 - 4 

3 

4- 6 

- 9 

2 

-2 

2 

-6 - 5 

4 

1 

- 8 

3 

-1 

5 

-5 - 6 



- 7 

5 

0 


-4 - 9 

U2T+ 17"*6 

- 6 

8 

1 

4- 5 

--3 - 16 

6 

-0"*016 

- 5 

25 

2 

4- 2 

-2 - 29 



- 4 

-f 33 

3 

9 

-1 - 68 


2D-^ 

- 3 

4- 10 



0 [-. 11] 

-13 

-0"-001 

- 2 

4- 5 


2/ -2D 

1 4- 91 

-12 

2 

- 1 

4- 3 

0 

- 0"-002 

2 4- 64 

-11 

2 

0 

-h 1 

1 

4 

3 - 127 

-10 

3 



2 

11 

4 - 7 

- 9 

6 

2D- 

27'-18°-2 

3 

76 

5 - 1 

- 8 

8 

— 5 

- 0"*004 

4 

4- 3 


- 7 

-' 13 

-6 

83 



2 4-7^-88“ 

- G 

22 

-7 

4 - 3 

2/-2D4-T-88* 

0 -0"-002 

- 5 

39 



4 

- 0"-004 

1 - 8 

- 4 

87 


2D4-? 



2 + 13 

- 3 

- 716 

-9 

4 - O^'OOl 

2^-2r 

» + 2r4-i7"-r) 

3 4 - 6 

- 2 

4 - 152 

-8 

-h 2 

6 

4 - 0"-065 

4 4 - 8 

1 

4- 74 

-7 

4- 2 



5 - 4 

0 

+ 39 

-6 

4- 2 


3/ -2D 

6 - 2 

1 

-H 13 

-5 

4- 3 

3 

-0"*003 


2 

4 - 10 

-4 

4 - 1 



/-7^-h88‘’ 

3 

4- 7 

-3 

4 


2Z-4D 

o 

o 

( 

1 

4 

4 - 5 

-2 

-f 4 

3 

4 - 0"-008 

-4 ^ 2 

5 

4 - 3 

- 1 

4 - 3 



-3 - 5 

6 

2 

0 

4- 2 


-4D 

- 2 - 13 



1 

+ I 

3 

- 1 - 0"-007 

-14- 8 

2D- 

l + T-60° 





0 + 2 

-1 

4-0" 002 




2F-2D 






3 

- 0"-002 

5r + 100° 







3 -0"-007 





2F- 

20 4-7^-90" 






5 

4- 0"*054 


2P-2D+; + ?'-90“ 
5 + 0"003 



ADDENDUM. FINAL VALUES 
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Final coefficients of sines of i {M — T) + arguments at top of each column. 


i 0 

2D -4/ + 149“ 

2D-Z 

1 +0"019 

- 5 + 0"-003 

- 6 - 0"*001 

2 + 8 


- 5 - 3 

3 + 3 

2D-24/+ 117“ 

-4 - 17 

4 + 2 

- 6 - 0"*002 

-3 + 5 

5 + 1 

±l 

-2 + 2 

+212“-7 

- 1 -- 0"-002 

2D~^-J/+ 149“ 

1 + 0"-002 


- 7 + 0"-001 

0 + 5 

Z+4/+ 212“ 

-<6 + 3 

1 + 60 

1 + 0"*007 

-5 + 38 

2 - 10 

2 - 1 

-4 - 3 

3 - 4 

4 - 2 

5 + 8 

->3 - 1 

5 - 1 

Z-4/+ 147“ 

2D-/^-24/+ 117“ 


- 2 + 0"-001 

- 6 - 0"-030 


-1 - 7 


1 -0"-002 


2Z-2D 

2 - 12 

Z + 24f+63“ 

4 - 0"‘002 

3 + 5 

2 -0"001 


4 + 2 

0 + 6 

2^-2D + 4f+2ir 

5 + 1 

^-2il/-63“ 

5 +0"-015 

34f + 96‘’ 

-2 +0"-001 

2^-2D + 24f + 63“ 

3 - 0"-002 

4 + 1 


6 -0"019 

3Z-2D + 24f + 63 

6 - 0"*001 
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ADDENDUM. FINAL VALUES 


Final coefficients of sines of i (J 

i 0 I ' 

1 +0''*069 -3 +0"-001 

2 - 13 -2 4- 4 

-1 + 8 
0 

0 -0"-209 1 ~ 8 

1+11 2 171 

2 + 8 

l+]b^ 

2J-W 3 +0"-002 

0 -0"-009 

2D 0 -.0"-021 

--5 + 0"-001 1 + 2 

- 4 + 2 2 + 58 

-3 + 3 3 - 1 

-.2 - 45 

-1 + 2 /-J+6° 

0 + 2 -1 - 0"’002 

1 0 + 21 , 

2 - 2 

+ 106'’ 

2D + J-6'’ 3 -0"-008 

0 - 0"-002 

/ + 2/-18° 
2 D - 7 + 7 “ 2 + 0"-002 

- 2 + 0"-020 

0 + 2 h + J-h" 

0 -0''*004 


T) + arguments at top of each column. 


2D-/ 

21 

-4 -0"-004 
-3 - 20 

2 -0"-013 

- 2 - 804 

2UJ-r 

-1 + 7 

2 +0"-002 

0 + 3 

2/- 2D 

2D-/ + J-100° 

2 -0"-187 

-3 +0"-004 

2Z-2D + J-6°‘7 

21)-/ + /^ 7° 

•2 -0"-190 

-2 +0"-007 

2^-2D + 2J-18' 

2D-/-/+ 6°73 

2 -0"'002 

-4 -0''*001 

-3 - 6 

2/-4D 

- 2 + 306 

2 +0"-009 

2D-/-y-80° 

2/-4D + J-7“ 

-1 +0"-005 

2 -0"-006 

2D-/~2./+18'’ 

3^-2D 

-2 +0"009 

2 -0"-007 

2D-/-2J^+107" 

3Z-2D + ^-7^ 

- 1 + 0"-002 

2 -0"-005 

2D + /5 

2i^-2D 

-2 -0''*003 

2 -0"*002 

2 - 1 

2F+l-2T> 

2D+Z-/+7‘’ 

2 -0"*002 

-2 +0"*001 

-2i^ + Z-2D 

4D-^ 

2 +0"-001 

O 

o 

o 

I 

1 


w - i - j+r 


-2 +0”003 




ERRATA. 


Page 10, line 7'| 

„ 13, „ 7 [ I'V --00706612- -00814813 rea<Z + -006624 <2+ •00826013. 

„ 14, „ loJ 

Page 10, line 8^ 

„ 13, ,, 8 [ l^’or - -008148 13+ -001210/3 mid f 008260 Jo- -001238 13. 

„ 14, „ lli 

Page 14, line 3 from foot, For - -02161 read - *00812. 

Page 15, line 7, For - *4753 read + *4458. 

Page 16, lines 4, 5, For round its centre of mass read round the centre of mass of the 
earth and moon. 

Page 19, equation (18), For R= read ii = real part of. 

Page 69, line 2, For dD-F + l + 2T-SQ read 3D - 1 + 2T -SQ. 

Page 75, line 11, For +722 read +726. 

„ „ „ -0"-007 „ -0"-003. 

„ M n -147 „ -148. 


Page 77, column headed <pi) = 1 + 3T - lOK, change the minus signs before w\ w'\ /i" to 
plus signs. 

Page 79, line 9, 1st column. For +86 read +83. 

„ „ 10, „ „ + 5 „ + 2. 

,, ,, 8 from foot, For 0-2F read 0 + 2F. 

Page 80, line 1, For 3h + 24:T read 3/i + 22T. 

„ ,, 4 from foot, 2nd column, For +0"-029 read +0"-027. 

,, ,, 3 ,, ,, Delete 0 + / +0"-002. 

„ „ 5 „ For -0"-0152 read -0"-0086. 

,, last line, For (p -0"-013 read 0 + 110^ +0"-007. 


cambbidoe: printed by john clay, m.a. at the university press. 








